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We review recent findings on spin glass models. Both the equilibrium 
properties and the dynamic properties are covered. We focus on progress 
in theoretical, in particular numerical, studies, while its relationship to 
real magnetic materials is also mentioned. 



The motivation that pulls researchers toward spin glasses is possibly not 
the potential future use of spin glass materials in "practical" applications. 
It is rather the expectation that there must be something very fundamental 
in systems with randomness and frustration. It seems that this expectation 
has been acquiring a firmer ground as the difficulty of the problem is appre- 
ciated more clearly. For example, a close relationship has been established 
between spin glass problems and a class of optimization problems known to 
be iVP-hard. It is now widely accepted that a good (i.e. polynomial) compu- 
tational algorithm for solving any one of iVP-hard problems most probably 
does not exist. Therefore, it is quite natural to expect that the Ising spin 
glass problem, as one of the problem in the class, would be very hard to 
solve, which indeed turns out to be the case in a vast number of numerical 
studies. None the less, it is often the case that only numerical studies can 
decide whether a certain hypothetical picture applies to a given instance. 
Consequently, the major part of what we present below is inevitably a de- 
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scription of the current status of spin glass studies, what is being said and 
done on the issue, rather than an account of estabhshed facts. In what 
follows, we put an emphasis on the theoretical and the numerical progress 
achieved in the field in the last decade. We mention older theoretical and 
experimental works only when it is necessary for readers to be able to follow 
current topics without reading through many other articles. Many of the 
important older works and experimental results, therefore, have been left 
out. For a review on these, the readers are referred to P, [2], 0, and 

The present paper is organized as follows. In section ^ we first give a 
brief overview of two well-known paradigms on spin glasses and summarize 
the predictions derived from them. Then, in the subsequent sections 12 
|31 and 01 we discuss equilibrium properties of Ising spin glass models. In 
particular, in section |21 a recent active debate, as to which paradigm is 
appropriate for realistic short-range spin glass models in three dimensions, 
is presented. Then, we proceed with the dynamical properties in section 
|5l to examine these paradigms on a different ground. In this section, an 
emphasis is put on aging phenomena while some other non-equilibrium 
properties are also discussed. Models with continuous degrees of freedom 
as well as the Potts spin glass models are mentioned in sectional The effects 
of weak disorder discussed in section provide us with a set of interesting 
issues quite different from those arising from strong disorder dealt with in 
the preceding sections. Several exact relations can be found and play an 
important role in shedding light on the issue under debate. Finally, we see 
in section IHl what results from the interplay between quantum fluctuations 
and randomness in spin glasses. 

1. Two Pictures 

The equilibrium properties of Ising spin-glass models in finite dimensions 
have been investigated mainly through developments of working hypotheses 
and numerical techniques. In principle, we can judge which working hypoth- 
esis is correct by numerical computations. The working hypotheses are often 
called "pictures" . Among various pictures, unarguably the most frequently 
mentioned are the replica-symmetry-breaking (RSB) and the droplet pic- 
ture. In fact, the dichotomy between these two has been the central issue 
in the spin glass research for more than a decade. However, it is rather 
difficult to perform a conclusive numerical test in the three dimensional 
case. One thing that renders the test so difficult is that the lower critical 
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dimension appears to be close to three. Another factor is, as mentioned in 
the introduction, that there is no polynomial algorithm for obtaining the 
ground state of a given instance in three (and higher) dimension. Therefore, 
the best we can do at present is to summarize on-going arguments, discrim- 
inating what has been established from what has not, and presenting the 
most important numerical evidences. We try to do these in what follows. 

It should also be pointed out that most of the "pictures" are asymptotic 
theories which are supposed to be correct only for large systems. In contrast 
to most regular systems the utility of an asymptotic theory for glassy sys- 
tems cannot be taken for granted. It was pointed out^ that in most real spin 
glass materials the typical length scale of an equilibrated domain within the 
experimental time scale may not exceed a small number of lattice spacings. 
If this is the case, considering the numerical evidence of the existence of 
large corrections to scaling, it can happen that the asymptotically correct 
theory fails to explain experimental results. The fact that many experiments 
on spin glasses, those at temperatures below the transition in particular, 
are not performed in equilibrium implies that an equilibrium theory may 
not be appropriate. Some attempts to construct a theory for a realistic time 
and length scales are reviewed in section |5| 

1.1. Mean- Field Picture 

We start with a brief description of the two pictures mentioned above in the 
present and the next subsection. We do not, however, intend to repeat the 
whole history of the development of each picture or discuss every technical 
detail. The reader may find a complete description of these pictures in a 
number of review articlesP^EEEGSl 

The RSB picture or the mean-field picture is based on Parisi's 
solutioJnUaSEl of the Sherrington-Kirkpatric (SK) modeP and its in- 
terpretation in terms of a multitude of thermodynamic states |16I17| . The 
SK model of N Ising spins. Si, is defined by the Hamiltonian 



where the summation is taken over all N(N — l)/2 pairs of spins with 
i j. The coupling constants Jij are quenched Gaussian random variables 
with zero mean and variance J/N where J is a positive constant of 0(1) 
that sets the energy scale. In jTH], the replica method was used in which 
n identical copies (replicas) of the system are introduced to perform the 
disorder average of the free energy: —(3F — [lnZ]av — hm„^o([2^"]av — l)/^i- 
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They found a solution of the mean-field equations using a replica symmetric 
Ansatz (see below). This solution has a phase transition and the transition 
point is located at the temperature Tc = J. Below this temperature, the 
order parameter 



has a non-zero value where indices a and (3 specify replicas. The replica 
symmetric Ansatz means that only solutions to the mean-field equations 
with q"^^ being independent of a and (3 (before taking the limit n — s- 0) are 
allowed. However, this solution did not describe the low temperature phase 
correctly since the entropy predicted by this solution became negative at 
low temperatures. 

In a different approach was taken for the SK model. They con- 
structed a set of equations (called the TAP equations) that expresses the 
equation of state for each bond-realization of the SK model in terms of 
the local magnetization, rui — (Si). They studied the eigen-modes of the 
stability matrix of stable solutions of the TAP equations and foimd that 
the spectrum of the eigen-values extends down to zero, suggesting that the 
solutions are only marginally stable in the thermodynamic limit. Further- 
more, the stability matrix in the replica space, i.e., the second derivatives 
of the free energy with respect to the replica order parameters, q"^, was 
examined. It turned out that the replica symmetric solution is stable only 
in the paramagnetic phase and becomes unstable below the phase boundary 
called the AT line in the T-H plane. Therefore, the SK solution is correct 
above the AT line whereas it is not below the AT line. 

To find the correct solution of the SK model below the AT line the 
replica symmetry had to be broken. The way the symmetry is to be broken 
is, however, highly non-trivial. A solution, called the Parisi solution, was 
proposec^E^EEI based on a novel metric structure of the replica-index 
space. It has now been proved with acceptable mathematical rigor that the 
solution is exact.l^Slln the Parisi solution, a metric structure is introduced 
in the replica-index space and it is assumed that the overlap q between two 
replicas depends only on the distance between the two in this index space. 
As a result, the order parameter is actually not a single number, in contrast 
to the EA order parameter qEA, but a function q{x) (0 < x < 1), where x 
is a suitable parameterization of the distance between two replica indices. 
A remarkable feature of Parisi's solution below the AT line was that q{x) 
varies continuously as a function of x. 




(1) 
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The order-parameter function q(x) first appeared to be a mathemati- 
cal artifact, but its physical meaning was later clarifiec^^by relating q{x) 
to the probability distribution P{q) of the overlap q between pure states, 
rather than the replicas. A pure state is defined as an extremal equilibrium 
distribution in the thermodynamic limit that cannot be expressed as a lin- 
ear combination of any other distributions. Denoting the inverse function 
of q{x) as x{q), it was found that 

x{q) = r dq'P{q'). (2) 

Here, the overlap distribution P{q) is formally (and somewhat symbolically) 
defined as 

P{q)=Y^P^Pp5{q~q^^). (3) 

Here the indices a and (3 specify pure states, not replicas, qap is the overlap 
between the two pure states a and /3, and Pa is the probability for the 
system to be in the pure state a. 

Equation |(2Jl may provide a way to investigate the structure of the space 
of the equilibrium distributions with the replica method. We can conclude, 
for instance, that more than one pure state exists and the overlap between 
two pure states varies depending on the particular pair of states if q{x) is 
not constant. If q{x) changes continuously in a finite range of x, as in Parisi's 
solution for the SK model, there must even be an infinite number of pure 
states. Simply stated, the RSB picture for systems in finite dimensions with 
short range interactions consists of the hypothesis that there are infinitely 
many pure states with varying overlaps and various predictions derived 
from the hypothesis. 

The definition Q of P{q), however, must be taken with caution since 
the existence of pure states as a well-defined thermodynamic limit is not 
clear for disordered systems. In fact, it was showiPE3that Pj{q) for a bond 
realization J depends on J even in the thermodynamic limit. As discussed 
below, this meanP^that a unique thermodynamic limit of Pj{q) with fixed 
bond reahzation J does not exist. Therefore, the order parameter, q{x), is 
not self-averaging for the SK model, in striking contrast to homogeneous 
systems. 

Because of this subtle nature of the order parameter, the exact implica- 
tion of the RSB picture for more realistic spin glass models in finite dimen- 
sions is somewhat ambiguous and has been a subject of debates.l^ESl One 
focus is the interpretation of the numerical result for finite systems. The 
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definition of q and P{q) often used in numerical simulations is based on the 
overlap between two spin configurations, S and S' , independently chosen 
from the equilibrium distribution of the same bond realization. Namely, 

1 ^ 

1=1 

Pj{q; i) = E E ^^(^; L)Pj{S'; L)5 {q - q{S, S')) 
S S' 

where Pj{S;L) is the normalized Boltzmann weight of the spin config- 
uration S* in a system of the linear size L with the bond configuration 
J. Throughout this article, unless stated otherwise, the angular bracket 
[• • • ] denotes the bond-configuration average with the distribution P{J) — 
Y[{ij) where the single-bond distribution Pi is independent of (ij). 

The following link-overlap and its distribution are also often discussed re- 
cently: 

qi{S, S') = ^Y. SrSjSlS!^ , P{qi;L) = [Pj{qi;L)] , (5) 
' (y) 

Pj{qi;L) = E E L)Pj{S'; L)5 {qi ~ qi{S, S')) , 

S S' 

where the summation concerning (ij) is over all the pairs of nearest neighbor 
spins, and Ni is the total number of the pairs in the system. 

The following is the list of features that follows from the RSB picture. 

(1) The overlap distribution, P{q;L), converges to a function P{q) that 
has a continuous part in the limit of i — > oo. In particular, P(0]L) 
converges to a finite value P{0) > 0. 

(2) Some quantities, such as Pj{q', L), are non-self-averaging, i.e., a unique 
thermodynamic limit does not exist for a fixed bond realization. 

(3) The bond-averaged link-overlap distribution P{qi] L) is not 5{qi) in the 
thermodynamic limit, i.e., its width converges to a finite value. 

(4) Global excitations with the energy cost of 0(1) exist. 

(5) A change in the boundary conditions generally affects spins located far 
away from the boundary. 

1.2. Droplet Picture 

Another well-known picturPEHElEaini is based on a scaling hypothesis 
on local excitations and produces markedly different conclusion for various 
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quantities, in particular, P{q). The basic assumption in this picture is that 
even below the transition temperature in zero magnetic field there are only 
two pure states that are mapped to each other by the total inversion of 
spins as in the homogeneous ferromagnets. Then, several scaling properties 
are assumed for compact excitations with varying size, which are called 
droplets. To be specific, a droplet of the scale I has a typical excitation 
free energy ci ^ Tl^, where 9 is called the droplet exponent. Since this 
exponent relates the energy scale to the length, similar to the stiffness 
exponent that relates the domain wall excitation energy to the system 
size, the simplest assumption is to identify 6 with the stiffness exponent 0s ■ 
Among various results derived from this picture, of particular impor- 
tance are the scaling forms of the correlation functions F^ I ^'^ l For example, 
it is predicted for models with a continuous bond distribution and with no 
magnetic field that the following asymptotic form should apply: 

[{S,S,f]-[(SmS,r]-T/TR' (6) 

near zero temperature, where R is the distance between site i and site j and 
the thermal average (• ■ ■ ) is taken in a single pure state (if there are two). 
It follows that the variance of the distribution function P{q; L) defined in 
(0J (more strictly, the variance of P{\q\;L)) decreases as the system size 
increases: 

(Ag)2^[(g2)]-[(|5|)]2«i-«^0. 

Therefore, when P{q) is defined as the thermodynamic limit of P{q;L), it 
consists of a pair of delta peaks at g = i^EA- 

For the models with discrete energy levels, such as the ± J model, the 
scaling form of the correlation function near T = must be modified. 
However, it is widely believed that the discretized nature of the energy 
level is not relevant at finite temperatures and the behavior is qualitatively 
the same as the one of the models with continuous energy levels. 

Below, we list the defining properties of the droplet picture together 
with some predictions derived from it. 

(1) The distribution of the excitation free energy of droplets of the scale /, 
P(e;), is continuous down to zero energy. 

(2) P{ei) has the typical energy scale Tl^ with T being an 0(1) constant. 
Specifically, P{ei) has the scaling form 
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(3) P{q) is self-averaging and it consists of a pair of delta peaks below the 
transition temperature. 

(4) Since P(0; L) is proportional to the excitation probability of a droplet 
that contains approximately one half of all the spins, for a system with 
a continuous bond distribution, P(0; L) has the following temperature 
and size dependence near T — 0; 

P(0; L) - [Prob (e^ < T)] - ^ cx TL'^. (7) 

EL 

(5) For systems with a continuous bond distribution the bond-averaged 
link-overlap distribution P{qi;L) consists of a single delta function in 
the thermodynamic limit. At low temperature, its width depends on 
the size and the temperature a^^ 

(Aqi)^ oc TL-^ i^ii = 2ds -2d~e) (8) 

where dg is the fractal dimension of the surface of droplets. 

(6) For systems with finite temperature spin glass transition, 9 is positive. 
Therefore, a global excitation costs an infinite energy. 

(7) A change in the boundary conditions does not affect spins located far 
away from the boundary, except for a potential simultaneous inversion 
of all spins. 

2. Equilibrium Properties of Two-Dimensional Ising Spin 
Glasses 

Two-dimensional spin glass models are easier to study than three- 
dimensional models from a technical point of view. However, there are 
still some open questions. Among others, the possibility of the existence 
of a phase transition at non-zero temperatures in the discrete distribution 
models is important. 

2.1. Zero- Temperature Transition? 

One of the arguments that supports the zero-temperature phase transi- 
tion scenario is based on the numerical estimates of the stiffness exponent 
that characterizes the size dependence of the domain- wall free-energy. The 
stiffness is simply the difference in the free energy caused by a particular 
change in the boundary conditions. The most frequent choice in numerical 
computations is the application of periodic versus anti-periodic boundary 
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condition. The stiffness exponent 9s is then determined by the asymptotic 
size-dependence of the stiffness, AF: 

AFcx 

According to [53! ; when the distribution of couphng constants is continu- 
ous and has a non-zero weight at J = 0, whether a phase transition takes 
place at a finite temperature or not is determined by the sign of the stiff- 
ness exponent, ff it is positive, the zero temperature phase is strongly or- 
dered, leading to a finite-temperature phase transition, whereas otherwise 
the order is fragile and infinitesimal thermal fluctuations would destroy it, 
leading to criticality at zero temperature. Numerical computations of the 
ground states of finite systems in two dimensions were carried 
for estimating the stiffness exponent at T = 0. For the Gaussian distri- 
bution of the bonds, it turned out that 9s ^ —0.3. At present, one of the 
most accurate and reliable estimates of the stiffness exponent is given inj3()| 
as 9s = —0.282(2) for the Gaussian bond distribution with periodic/anti- 
periodic boundary condition. Due to the negative value of the stiffness ex- 
ponent, it is widely believed that there is no finite temperature transition 
in this model and that the system is critical at zero-temperature. 

The stiffness exponent of ± J model, on the other hand, is much closer 
to the marginal value, 0. It was estimatedj^ with the assumption of the 
power-low dependence of the stiffness on the system size, as 

9s = -0.060(4), 

while the possibility of 6*5 = was not ruled out. The possibility of the 
stiffness exponent being zero was strongly suggested by another ground 
state computation!^ 

Recently, it was pointed out that the stiffness exponent have to be inter- 
preted differently for ±J models. Amoruso et aL^^ carried out a renormal- 
ization calculation using the Migdal-Kadanoff method. They found that the 
stiffness exponent, calculated within this approximation, is zero in any di- 
mension lower than the lower critical dimension. This finding suggests that 
the stiffness exponent being zero does not necessarily imply that the system 
is exactly at the lower critical dimension. This type of dimensional depen- 
dence of the stiffness exponent was observed only for the class of the bond 
distribution for which the discretized nature of the coupling constant is not 
smeared out by the renormalization. Since the ± J bond distribution falls 
into this class, the above-mentioned numerical estimates, ^ 0, may only 
mean that the lower-critical dimension is two or above. For the ± J model. 
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however, even the possibiUty of a finite temperature phase transition was 
suggested|22i as we discuss in greater detail further below. 

2.2. Droplet Argument for Gaussian- Coupling Models 

For two dimensional systems with a symmetric {P{ J) ~ P{—J)) and contin- 
uous bond distribution, it was arguecPSl that there is only one independent 
critical exponent. In other words, all critical indices are related to the stiff- 
ness exponent via scaling laws. We consider the following finite size scaling 
form of the singular part of the free energy Fs, 

^{T,H,L) = = -[logZ] ^ (j){TLy,HLy^). 

The condition that the total magnetization [M] = lim/f^o dF^/dH is pro- 
portional to L"^/"^ at T = relates y to yh via 

d 

yh = y+^- 

Then, by differentiating $ with respect to T and H, we can express any 
critical index by y. For example, the magnetization per spin depends on 
the magnetic field as to = L~ M ^— L"^ ^ FL'^^h ^ which means that the 
exponent 8 appearing in the scaling of the magnetization with the magnetic 
field is (5 = 1 — lyjd. The non-linear susceptibility X2 = L~'^d^M jdFl^ 
at H — depends on the temperature as X2 ^ T~'^^ with 72 = 3 + 
^. Since the spin glass susceptibility XSG at iJ = is related to X2 as 
XSG = L-'^J2i,j{{Si^jy] = T^X2, the corresponding exponent is 730 = 
d/y. Similarly, the specific heat exponent can be expressed as a = —d/y. 
As for the exponent ?7sg: that characterizes the asymptotic form of the 
two-point spin glass correlation function [(5*^5^)^] ~ R'l~'^'^^^^ , we have 
2 — ''7SG = 2/7sGj which yields jysG = 2 — d. This is consistent with the fact 
that the correlation function does not decay at zero-temperature. 

2.3. Droplets in Gaussian- Coupling Models: Numerics 

As mentioned above, one naturally expect^" that the exponent —y is iden- 
tical to the stiffness exponent 9s since both the exponents relate the en- 
ergy scale to the length scale. However, the estimates of various critical 
exponentPEsEnEHEH 

based on the system size up to L 50 seem to sat- 
isfy the scaling relation with y ~ 0.5 rather than y = —9s ~ 0.3. 

In order to check if y coincides with —9, a direct numerical estimate 
of the droplet excitation exponent 9 was carried out using a heuristic op- 
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timization procedur^^El applied to the EA model with a Gaussian bond- 
distribution in two dimensions. For each realization of the model, the ground 
state spin-configuration was computed with free boundary condition. Then, 
the spins on the boundary were fixed as they were in the original ground 
state whereas the spin at the center was fixed in the opposite direction. 
These constraints lead to a new ground state that is identical to the old 
one on the boundary but differs from it in the vicinity of the central spin by 
a droplet of flipped spins around the center. The droplet is typically as large 
as the system itself. The system-size dependence of the droplet volume V 
and the droplet excitation energy E could be described well by the scaling 
law 

l^(L)cxL'^°, and E{L) oc {V{L)f/'^'' 

with 

dr, 1.80(2), and - 6* = 0.47(5) 

for the range of the system size 5 < L < 49. In particular, the value of —9 
agreed with most of the previous estimates of y. This can be understood 
also within the droplet theory because all scaling forms derived by droplet 
arguments are identical to what one can get via the ordinary finite size 
scaling by identifying 9 and —y. 

However, a recent computatiorPl demonstrates the presence of a cross- 
over in the droplet excitation energy. They performed essentially the same 
calculations as the one described above for larger systems (up to L = 160). 
and argued that there may be a correction term due to the self-interaction 
of the surface of the droplets and that the droplet excitation energy has the 
form 

AE{L) - AL^ + BL-'^ 

with Lu > \9u\. The numerical data could be well fitted by 

-9d = 0.29, and lu = 0.97(5). 

as is shown in FigQ^a). The fact that the estimate of the droplet ex- 
ponent now becomes close to the stiffness exponent indicates the validity 
of the simplest scenario. However, the result on the fractal dimension of 
droplets seems to suggest the contorary; their result appears to confirm the 
previous result do ^ 1.8 (Fig^b)). While this may be another transient 
behavior, the problem of compactness of the droplets remains open in the 
two-dimensional system. 
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(a) 



(b) 
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Fig. 1. (a) The droplet energy as a function of the system size. The two dashed straight 
lines represent the algebraic dependences L"'' *^ for smaller systems and L~^-^^ for larger 
ones, respectively. The solid curve is a fitting function (see text), (b) The volume and 
surface of the droplets. (From f42| .') 

The same model was studied in 03] and then in 0^ at zero-temperature 
using a different approach. For each reaHzation of the model, the ground 
state for periodic boundary condition was compared with the ground state 
for another set of boundary conditions, such as anti-periodic and random 
boundary conditions. (In the latter, half of the bonds across the bound- 
ary are chosen at random and inverted while the other half are kept un- 
changed.) The authors of 03] argued that in the RSB picture the change 
in the boundary condition would induce an excitation whose surface-to- 
volume ratio should not diminish as the system size increases, whereas in 
the droplet picture it should decrease down to zero. Therefore, we can dif- 
ferentiate between the two pictures by applying a boundary condition that 
induces a domain wall in the system and focusing on a small box located 
at the center of the system. If the probability of having the domain wall 
crossing this part is finite even in the thermodynamic limit, it would indi- 
cate the validity of the RSB picture. On the other hand, it should decrease 
to zero if the droplet picture is valid. 

They performed a calculation of exact ground states of many instances 
of the EA model in two dimensions with the Gaussian bond distribution. 
The system size explored was up to i = 30. They considered a box of the 
linear size of Lbox = 2 located at the center of the system, and measured 
the distribution of the overlap qbox on this box: 
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where a and f3 are the indices specifying real rephcas. If the boundary 
passes the box with a finite probabiUty, the distribution of should have 
a finite weight around q^^^ = 0. This weight turned out to be decreasing 
with the power of —0.7 as a function of the system size. Based on this result, 
they concluded that the zero-temperature structure of the two-dimensional 
EA model is trivial, i.e., only two pure states exist. (See Fig0]) 

2.4. Finite- Temperature Transition? 

There are some numerical results^ of Monte Carlo simulation on the ±J 
model that suggest the existence of a phase transition at a finite tempera- 
ture. Specifically, the overlap distribution P{q; L) and the binder parameter 
g were computed. The estimates of the Binder parameter 



a finite transition temperature ~ 0.24J was assumed in stead of zero 
transition temperature. Their estimate for the critical indices of this finite 
temperature transition are v ^ 1.8 and rj ^ 0.2 The overlap distribution 
P{q; L) could be nicely scaled with the same rj but with a slightly higher 
temperature Tc ~ 0.29 J. Similar results were obtainea^'^ for an asymmetric 
bond distribution. 

Numerical evidences suggest that the stiffness exponent of the two- 
dimensional ± J model is non-positive, as discussed in subsection l2.1l which 
indicates that the lower critical dimension of the ± J model is equal to or 
greater than two. Now, if a spin-glass transition takes place at a finite tem- 
perature as suggested, it means that the lower critical dimension is not 
greater than two. Therefore the only scenario consistent with all the avail- 
able numerical and analytical results is that the lower-critical dimension is 
exactly two and that the stiffness exponent is exactly zero in two dimen- 
sions. More evidences, however, appear to be necessary to settle this issue 
beyond reasonable doubts. 

3. Equilibrium Properties of Three-Dimensional Models 

The problem of spin glasses in three dimensions is the central topic of 
the field. Whereas the RSB nature of the low-temperature phase in four 
dimensions is much less controversial, the nature of the low-temperature 
phase in three dimensions still remains the subject of an active debate. 




(9) 



as a function of the temperature for various system sizes appeared to have 
a common crossing point. The finite-size-scaling plot seems better when 
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Below we present a number of numerical results. Because of the severe 
technical limitations for three dimensional cases, all the results obtained 
are for small systems, typically up to the linear size of about 10 lattice 
spacings or slightly more. Therefore, many important issues are left open, 
in particular the question as to what picture yields the correct description 
of the low-temperature phase of the EA model, which is the main subject 
of subsection \'A . 31 through subsection 13. 81 

3.1. Finite Temperature Transition? 

Even the very existence of a phase transition at finite temperature was 
not easy to establish. An evidence of the existence was obtained through 
the calculation of the domain wall energy at zero-temperaturCjI^Sl in which 
a positive estimate for the stiffness exponent was obtained. However, the 
system size was rather limited {L < 8) and the estimated value of the 
stiffness exponent was small {6 ^ 0.2). Therefore, this finding about the 
stiffness exponent alone was not sufficient to establish the existence of a 
finite-temperature phase transition. Finite temperature approaches could 
not settle the issue, either. For example, while Monte Carlo simulationPsEI 
strongly suggested the existence of a transition at a finite temperature, they 
could not rule out the possibility of zero-temperature singularity with an 
exponentially diverging correlation length; it was suggested'** that all exist- 
ing data at that time (namely 1994) were consistent with both hypotheses: 
Tc > and Tc = 0. In particular, the simulation results of for the spin- 
glass susceptibility could be fitted by a functional form consistent with a 
zero-temperature singularity, XSG = l + A/ (T — Tc)'^, just as well as the one 
consistent with a finite-temperature transition xsG = 

A {e^B/rr -i)+c. 

However, another set of Monte-Carlo resultPl showed beyond reason- 
able doubts the existence of a finite temperature transition in the ± J model 
in three dimensions. Specifically, it demonstrated that the Binder parameter 
curves for different system sizes cross at the same point (FiglSJ near 

T - 1.1, and ^^g ^ 0.75. 

With the assumption of the algebraic singularity at the critical point, the 
critical temperature and the critical indices were estimated as 

rc== 1.11(4), J/ = 1.7(3), and 77 = -0.35(5). 

These estimates were consistent with previous one^^ and confirmed by 
other simulations as presented in Tabled 
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Fig. 2. The Binder parameter vs temperature for the it J model in three dimensions. 
(From 49 .) 

A Monte Carlo simulation for larger system sizeJ^ clarified the issue 
further. The results of the correlation length and the spin-glass susceptibil- 
ity were extrapolated to the infinite system-size limit via finite size scaling, 
and were well fitted by a curve representing an algebraic divergence at a 
finite temperature. While the existence of a finite temperature phase tran- 
sition could be concluded from this result, the nature of the singularities 
at the critical point could not unambiguously be settled. (The singularity 
could be an essential rather than an algebraic singularity.) In the next sub- 
section, we present other numerical results concerning the nature of the 
critical point. 

3.2. Universality Class 

Early Monte Carlo simulations such as 0^1 indicated a scenario that the sys- 
tem stays critical in the whole low temperature region. Namely, these data 
could be explained by a line of critical points terminating at Tc = 1.2, simi- 
lar to the Berezinskii-Kosterlitz-Thouless transition in the two-dimensional 
XY ferromagnet. The conclusion of [2111 on this issue was similar, while it 
was suggested that a finite-temperature transition with a algebraic diver- 
gence is most likely. With the assumption of the algebraic divergence, the 
critical parameters were estimated to be consistent with previous ones men- 
tioned above. However, the data were also consistent with an exponential 
divergence as in the BKT transition. 
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16 N. Kawashima and H. Rieger 

A large-scale Monte Carlo simulatior^clarified the issue of the nature of 
the phase transition. Using special purpose machines and parallel tempering 
methocP^EH the authors succeeded to equilibrate the ±J model of sizes up 
to L = 20 down to a temperature low enough to cover a sufficiently large 
region around the transition point. They obtained various results consistent 
with an ordinary second order phase transition. For example, they observed 
a clear crossing in the effective correlation length-"* £,{L) divided by the 
system size L. (They defined the effective correlation length as ^^{L) = 
k^{Cq{Q)/Cq{km) — 1), whcrc Cq{k) is the static structure factor, km is 
the smallest non-zero wave number compatible to the boundary condition, 
and fc2 = 4(sin2(fc^/2) -I- s\T?{ky /2) + sin2(fcV2)).) 
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Fig. 3. The effective correlation length divided by the system size. The left panel is 
for the ±J model in three dimensions whereas the right panel is for the XY model in 
two dimensions. (From 1511 .) 

The scaling dimension of the effective correlation length is zero; the 
curves for various system sizes should intersects at the same point similar 
to the Binder parameter. Indeed, it shows a common intersection point as 
in the right panel of Fig|31 This is in marked contrast to what we typically 
see in a BKT type transition, which is shown in the right panel of Fig|3| 
for comparison. From this result one concludes that the phase transition is 
an ordinary second-order one. For the critical temperature and the indices, 
again, see Tabled 
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The results of provided also indications in favor of an RSB nature 
of the system below the critical temperature. Namely, the estimates of a 
quantity called the G-parameter seems to have a finite value even below 
the critical point. The G-parameter is definecP^as 

The numerator is finite only if the sample-to-sample dependence of P{q) 
does not vanish, while the denominator can be finite even if there is no 
sample dependence. Therefore, a observed non- vanishing value of G implies 
that P{q) is not self-averaging. 

However, there are a number of evidences by other groups that can be 
interpreted otherwise, namely, in favor of the absence of an RSB nature of 
the low temperature phase, as we discuss below. 

3.3. Low- Temperature Phase of the dzJ Model 

There is a very active debate on the nature of the low-temperature phase of 
three dimensional spin glass models. In the following few subsections, we re- 
view various theories, arguments and numerical calculations that were made 
or done with the ultimate aim to clarify the issue. But none of the pictures, 
RSB, droplet, or others, could be estabhshed so far. Nonetheless, the debate 
itself is interesting and each picture is worth being scrutinized with detailed 
numerical calculations. While the common belief is that the ± J model and 
the Gaussian-coupling model show essentially the same physics at finite 
temperature, we present them separately since many of the computations 
are done at zero temperature where the two may differ. The ±J model is 
discussed in the present subsection and the Gaussian-coupling model in the 
next. 

A finite temperature simulation on the ± J model was performeci^ using 
the replica exchange Monte Carlo simulatiorl^ for ± J model reaching down 
to r = 0.2 for the system size up to i = 10. The integrated overlap 
probability distribution 

.1/2 

x{l/2)= / dqP{q) 
Jq 

was computed as a function of the temperature and the system size. Inspired 
by EBj, the scaling — x{l/2) being proportional to for small L and to 
for larger L — was assumed. In other words, 

x{l/2) = TL-^fiTL^) 



Table 1. The estimates of tlie critical temperature and the exponents for the three dimensional spin glass models. The entries 
are categorized into three groups according to the anisotropy; easy-axis (top), isotropic (middle), and easy-plane (bottom). The 
entries with "*" are not quoted in the original paper but estimated through the scaling relations by the present authors. The 
entries with are for the chiral glass transition. . 
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was assumed. The fitting to the numerical results worked nicely with 9 = 
yielding 

A = 0.9(1), 

whereas a fitting with 6 = 0.2, the droplet prediction, turned out to be sig- 
nificantly worse. The estimate of A is considerably smaller than a preceding 
estimatJ^ while consistent with the value A = 0.72(12) in 80 . 

Zero-temperature calculations offer a powerful alternative to the finite 
temperature approach with Monte Carlo simulations, because one might 
expect that the differences between the two scenarios should be more promi- 
nent at lower temperature, and much better methods are available for sys- 
tems at zero-temperature than for a very low but finite temperature. For 
three dimensional models, most methods for solving zero-temperature prob- 
lems are based on heuristic optimization since no good exact method is 
available due to the NP-hardness of the problem. A computation of the 
ground states of the ± J Ising model up to L = 14 was donJSUHll with a 
heuristic algorithm called the cluster-exact approximation method.!^ Later 
the computation was redone ,'^^'^41 in order to fix the problem of the biased 
sampling!^ It was found that the width of the overlap distribution, P{q), 
decreases as the system becomes larger, indicating the triviality of P{q)- 
The validity of the ultra- metric relation, qi2 = (723 < 931, was also examined, 
where qij^s are the overlaps among three randomly chosen ground states. 
The numerical results indicated that the ultra-metric relation holds for a 
typical triplet of ground states with relatively large mutual distances. How- 
ever, it was also found that the contribution from these triplets to P{q), 
which constitutes the "non-trivial" (i.e., continuous) component of P{q), 
decreases as the system becomes larger. In fact, the integrated weight of 
P{q) systematically decreased toward zero as 

x{q) EE / dqP{q) cx 

J-q 

with A = 1.25(5) for q = 0.5. This indicated a trivial structure of P{q) at 
zero temperature. As for the exponent A, a different estimate was obtained^ 
as mentioned below. 

The authors of |2B] considered P{q) for the ±J model at T = 0, and 
pointed out that the ground state space may be dominated by a single valley 
or a single pure state no matter which picture is valid. They argued that 
even if there are multiple valleys in the energy landscape of the discrete 
energy model, the distribution of the overlap P{ci) may still be trivial, 
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indicating that it is impossible to discriminate between the two scenarios 
by a computation of P(ci) such as the one mentioned above. 

The underlying assumption of the argument is the many- valley structure 
in the phase space and the sponge-like structure of low-lying excitations in 
the real space. (The latter is further discussed below in subsection 13 . fil and 
subsection 13.71 ) A droplet typically has a finite volume and when flipped 
it takes a state in a valley to another state in the same valley. In contrast, 
a typical sponge-like cluster is supposed to occupy a finite fraction of the 
whole system and flipping it generally causes a transition from one valley 
to another. 

This argument is based on an estimate of the entropy by counting zero- 
energy excitations from a particular ground state. In general, the continuous 
part of P{q), if any, is caused by excitations of various scales whose excita- 
tion free energies are smaller than the temperature. In particular, at zero 
temperature, it is caused by zero-energy "excitations" . They argued that 
P{q) has a trivial structure even if there are multiple valleys. To see this, it 
suffices to consider the case of two valleys. Each valley contributes to P{q) 
according to its weight that is the number of distinct spin configurations in 
the valley. This number is roughly the same as the number of zero-energy 
droplets. We can choose two representative configurations, one from each 
valley, that can be transformed to each other by flipping a sponge-like clus- 
ter. Then, the zero-energy droplet excitations in one configuration differ 
from those in the other configuration only on (or near) the boundary of 
the sponge-like cluster. Therefore, if A is the area of the surface of the 
sponge-like cluster, the droplets differ only at A positions. This difference 
results in the entropy difference of ©(A^/^). (The power 1/2 comes from 
the assumption that the number of the droplet-like excitation located on 
the surface is a random variable.) It follows that the contribution (to the 
continuous part of P{q)) from one valley differs from that from the other 
typically by a factor e='=('^'^'^^^)^"^^^^. We now see that almost certainly P{q) 
is dominated by a contribution from a single valley in the thermodynamic 
limit where A ~^ oo, leading to a trivial structure in P{q). It is clear that 
the presence of more than two valleys does not affect the result as long as 
the number of valleys does not grow too fast as the system size increases. 
Based on this result, they argued that a trivial P{q) for the ± J model in 
three dimensions at T = does not necessarily indicate the absence of 
the RSB. 

However, it is not known how the number of valleys depends on the 
system size. It is not too unrealistic to assume that the number of valleys 
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grows faster than or proportional to the number of spins, N. If so, a typ- 
ical minimum entropy difference between two valleys would be less than 
A^^-^/N < 7V^^/^ 0. Therefore, it is not clear if this argument really in- 
validates the numerical evidence in |84| for triviality of the energy landscape 
structure. 

3.4. Low-temperature Phase of the Gaussian- Coupling 
Model 

The ground state of models with the Gaussian bond distribution, is unique 
up to a trivial degeneracy due to the Z2 symmetry. However, it is still pos- 
sible to extract useful information about low-lying excitations from zero- 
temperature computations. The authors of ,86^ attempted to use the same 
strategy that they used for two dimensional models with the Gaussian bond 
distributioi^ (see subsection 12 .3(1 to discriminate between the two scenar- 
ios in three dimensions. They estimated the probability of a domain- wall 
passing through a small imaginary box placed inside the system. A clear 
decreasing behavior as a function of the system size was found, in favor of 
the droplet picture. However, the amount of the total decrease that they 
could observe by changing the system size was only of a factor of 1.3 or 
1.4, due to a severe system size limitation. (In the case of two dimensions, 
the same quantity varies by almost an order of magnitude as can be seen 
in FigH) 




2 4 6 810 20 3 4 5 6 7 8 910 

L L 



Fig. 4. The probability of a domain wall passing through the box. The left panel is 
for the EA model with the Gaussian bond distribution in two dimensions at T = and 
H = 0, whereas the right panel is for the same model in three dimensions. (From 1861 .*) 
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In addition to the response to the boundary condition, the effect of a 
weak perturbation apphed to the whole system was also examined?^ The 
perturbation was chosen in favor of the excited states against the ground 
state. For each bond realization, they first obtained the ground state, which 
we denote {5'^^°^}. Then to each coupling constant Jij the following pertur- 
bation was added 



where AJ = e/iVbond was chosen to be of 0{L-'^). With this perturbation 
the ground state energy is increased, exactly, by e, whereas excited state en- 
ergies are changed according to its similarity to the ground state. Therefore 
by adjusting the amplitude of the perturbation and examining the overlap 
between the original ground state and the ground state of the perturbed 
system, one may obtain some information about the excitations. 

The RSB picture predicts the existence of many spin states that are 
local minima of the Hamiltonian. These states differ from the ground state 
by only an infinitesimal amount of energy per spin but differ by a macro- 
scopic number of spins such that the average Hamming distance between 
them (normalized by the system size) converges to some finite value in the 
large system size limit. This means for the bulk perturbation considered 
in j87j that the overlap between the two such minima, the one with the 
perturbation and the other without, should be truly smaller than unity in 
the thermodynamic limit, no matter how small the perturbation may be. 
Therefore in the RSB picture, 1 — [q\ should be finite where q is the overlap 
between the two minima. On the other hand, in the droplet picture, 1 — [q] 
should depend on the system size by a power-law being characterized by 
the droplet energy exponent 9 and an exponent related to the geometrical 
properties of the excitations. To be specific, since 1— [q] is dominated by the 
contribution from droplets whose linear size is comparable to the system 
size L, it is roughly equal to the volume fraction of a typical droplet of size 
L multiplied by the probability of such a droplet being excited. Namely, 



where is the energy gain by the droplet excitation due to the perturba- 
tion and is the droplet excitation energy without the perturbation. If the 
droplets are compact, as assumed in jHZ], only the second factor matters 
since = d. For A^,, the authors of [HT] assumed that e ~ 0(1) when inte- 
grated over the whole system, which is equivalent to setting A^ cx 




(10) 
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where ds is the fractal dimension of the boundary of droplets. Based on 
this, the following scaling form for the overlap was proposed j87| : 

v{e,L)^l~[q]^V{e/L''~''^-'). (11) 

For the link overlap qi, a similar argument yields 

s(e, L) = 1 - [qi] - L-('*-'^=)5(e/L'^-'*=-''). (12) 

In the system size range that could be explored, (namely, L < 8), the 
numerical results for [q] and [qi] shows that v{e,L) and s{e,L) obeys the 
power law and decreases to zero as the system size increases, contradicting 
the RSB picture. However, fitting the numerical results to the forms (|11(> 
and ifT^ yielded 

d-4 = 0.42(2), and d - 4 + 61 = 0.44(2). 

This implies 

e = 0.02(3). 

The difference between the value and the stiffness exponent 0.2) ob- 
tained from the domain wall calculation is statistically significant, i.e., the 
result contradicts not only to the RSB picture but also to the droplet pic- 
ture unless this seeming contradiction is caused by a large correction to 
scaling. 

Concerning the possible source of the correction, two argumentJ^^IllEll 
were presented to explain the inconsistency observed in numerical results 
within the framework of the droplet argument. It was pointed oulP^that if 
one assumes a "clean" scaling with no correction for the droplet excitation 
energy, the contribution from small droplets gives rise to a correction to 
scaling in quantities such as the spin-glass susceptibility and the magnetic 
susceptibility. The droplet excitation energy computed with the condition 
used in |40I41| and the link overlap computed in ^Wl\ may have been affected 
by such a correction to scaling. Another possible source of a correction to 
scaling is the interaction between domain walls. It was arguecP^EH that 
the energy of a domain wall may be increased by the presence of another 
domain wall and that the energy shift due to this interaction may have 
the form l^'^ where / is the distance between two domain walls. Similarly, 
a "self-interaction" of the domain wall may give rise to a correction to 
scaling in the droplet excitations a,s E = Al^ + Bl~'^ . In this case I stands 
for the size of the droplets. As the correction-to-scaling exponent due to 
this mechanism, the authors of [HH! quoted the value uj ~ 0.13(2^^. 
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In response to some reports (see the following paragraphs in the present 
section) contradicting to their conclusion, the authors of 91 further pur- 
sued the ground state nature of the three dimensional system with the 
Gaussian bond-distribution, along the same line as their own preceding 
calculation!S3 This time, however, larger systems {L = 12) were dealt with 
the branch-and-cut algorithm!^ This algorithm guarantees that the states 
found are the true ground states. In addition, a greater care was taken for 
various possibilities of fitting functions and different sources of corrections 
to scaling. Specifically, they considered a few different fitting functions for 
size dependent quantities such as the link overlap and the box overlap. Fit- 
ting functions consistent with the droplet, the TNT, and the RSB scenario, 
respectively, were considered. (For the TNT scenario, see subsection l3 . 61 and 
subsection l3.7l below.) It was found that the size dependence of the surface- 
to-volume ratio, i.e., 1 — q\, could be explained by any one of three pictures, 
and also that the size dependence of the box overlap could be explained by 
any one of three pictures. The estimates of the exponent /i = d — ds + 9 
turned out to depend on the boundary condition. With the free boundary 
condition, the estimates are d — ds = 0.44(3), /i = 0.63(3), and 9 — 0.19(6), 
whereas for the periodic boundary condition, they are d — ds — 0.43(2), 
fi = 0.42(3), and 9 = —0.01(3). The latter set of values are consistent with 
their preceding estimate whereas the former are not. 

Another set of numerical results, presented in |93l94l95j seems to con- 
tradict to the results of |87l86j presented above. As for the calculation based 
on the small imaginary box placed at the center of the system, it was ar- 
gued93l that the probability of the state inside the box being affected by 
the change in the boundary condition should obey the scaling 

^changc(-^box, -^) ^ g{Lhox/L), (13) 

if the droplet argument is correct. However, the numerical results in |93) 
for Lbox = 2, 3, 4 and L = 12 did not fit in this scaling. 

In I93j ■ the scaling property of the domain wall induced by the anti- 
periodic boundary condition in the x-direction was also examined. In par- 
ticular, the probability of the domain wall not intersecting a plane perpen- 
dicular to the X axis was measured. It was found that the non-intersecting 
probability approaches zero as the system size increases: 

PiC'The domain wall does not intersect the plane.") cx L~'^ , 

with J — 1.5 — 2.0. Based on this observation they suggested that the 
domain wall is space filling, i.e., d — ds- 
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In |94l95j ■ zero-temperature calculations were performed with the con- 
ditions analogous to the ones in and [HZ]- In IHlli the effect of the anti- 
periodic boundary conditions, which were imposed in the x direction, was 
compared with the periodic one. The links perpendicular to the yz plane 
were treated separately from those parallel to the yz plane. The result of 
the link overlap for the perpendicular links, qp, and that for the transverse 
ones, qx, are shown in FigEfa). Both kinds of the link overlap can be fitted 



(a) (b) 




0.05 0.1 015 0.2 0.25 '/ O-l (1.2 (I.IH 



Fig. 5. (a) The link-overlap for the perpendicular links, qp, and that for the transverse 
ones, q^, between the two ground states with and without a twist in the boundary, 
and (b) the link-overlap (subtracted from unity) between the two ground states with 
and without a bulk perturbation. (The lower set of data is for the average link-overlap, 
whereas the upper one is for the link-overlap restricted to those pairs of states which 
have zero mutual overlaps. For each set of data, the upper fitting curve is a second order 
polynomial in 1/L and the lower one is a fractional power in 1/L with an additional 
constant term.) For both (a) and (b), the model is the three-dimensional EA model with 
the Gaussian bond distribution at zero temperature. (From (a) '94' and (b) (^D.) 

well by second order polynomials in 1/L as shown in the figure. From the 
zero-th order term, they concluded that 

lim l-qi = 0.245(15), 

L — 'oo 

in contrast to the droplet prediction 1 — qi — s- 0. 

In (95| , the effect of the bulk perturbation at zero temperature was stud- 
ied. The perturbation was the same as the one considered in [HZ]. The link 
overlap, qi{q), was computed as a function of the bulk overlap, q. In addi- 
tion to the average link-overlap, qi{L), attention was paid to the quantity 
A{L) = l — qi{q = 0), where qi{q = 0) is the link overlap between two states 
that have zero overlap. In the droplet picture, both quantities should con- 
verge to as L goes to infinity. The numerical results are shown in FigE^b). 
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It was found that in both cases of the boundary and the bulk perturbations, 
the size dependence of the surface-to- volume ratio of the domain wall, i.e., 
1 — qi can be fitted by a polynomial in 1/L with a finite constant term, 
consistent with the RSB picture. 

Concerning the results presented in [HSI , it was pointed oulJ^Si that the 
probability Pchango is subject to a strong correction to scaling and that, 
taking such a correction term into account, the data presented in can 
be fitted reasonably well with the scaling form It was also pointed 

out that even the two dimensional data, which is believed to obey \V6\. 
could not be fitted by this scaling form without a similar correction term. 
As for the non-intersecting probability of the domain wall with the plane, 
it was simply pointed out that the observation in does not necessarily 
imply d = ds because the observed behavior can be caused by domain walls 
that are rough but not space filling. In addition, as already mentioned 
above, their latest calculatiorPD suggests that we have not yet accumulated 
numerical evidences sufficient to decide between the possible scenarios. 

Before concluding this subsection, let us mention what is known about 
the geometrical nature of droplets. For this issue, a zero-temperature calcu- 
lation was performed^ Using an efficient heuristic algorithm, the authors 
of 90 obtained ground states of three dimensional systems of the size L = 6 
and L = 10. For each sample, they first obtained the ground state. Then, a 
center of the excitation is chosen randomly. They searched for the droplets 
of various sizes including the central spin. In their calculation, a droplet was 
defined as the cluster of spins of the smallest excitation energy among those 
which have a given volume. In spite of the limited system size, they suc- 
ceeded in obtaining system-size independent numerical results. They found 
that the linear scale of the droplets is described by 

R{v)(xv'' (6-0.5), 

whereas the volume dependence of the excitation energy obeys 

E{v)(xv'' (a- -0.06). 

The first relation means that the droplets are roughly two dimensional 
fractal objects rather than compact ones. The second relation is rather 
surprising, since at first glance it might seem to contradict the existence 
of the finite-temperature phase transition. They claimed that this does 
not necessarily mean the absence of the phase transition because the basic 
assumption of the droplet argument may not be valid. Another possible 
source of the discrepancy may be the definition of the droplet with the very 
restrictive fixed- volume constraint. 
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3.5. Effect of Magnetic Fields 

Whereas the RSB picture predicts the persistence of the spin glass phase 
against a smah but finite magnetic field, the droplet-picture predicts the 
contrary — the absence of the spin glass phase at any finite magnetic field. 
Therefore, one of the possible ways to distinguish between the two scenarios 
may be found in the behavior of short-range models in a finite magnetic 
fields. However, it is rather hard to perform this task near the critical 
temperature (of the zero-field case). This difficulty is not only because the 
dynamics near the critical point is slow but also because the critical value 
of the magnetic field predicted by the mean-field type argument is small. 
To avoid at least the latter difficulty, a zero-temperature calculation was 
performed!^ Using a heuristic optimization algorithm, the ground state of 
the three dimensional EA model was computed with the Gaussian bond 
distribution up to L = 12 with various uniform magnetic fields. 

In addition to the magnetic field, also some constraints were imposed 
at the same time to probe low-lying excitations. For example, a spin may 
be forced to be opposite to its natural direction in the original ground 
state with no constraint. Then, the lowest-energy excitation among those 
which include the chosen spin is excited. In the RSB picture, where there 
are infinite number of states different from each other by a finite fraction 
of spins within an 0(1) window of excitation energy, such a excitation 
should typically contain a finite fraction of all spins. (See also subsection 
13. 61 and subsect ion l3 . 7l below . ) Therefore, the volume of excited cluster tends 
to infinity with increasing system size. On the other hand, in the droplet 
scenario, the influence of the forced spin cannot reach beyond the distance 
determined by the magnetic field H as 

liH) OC i/-2/(rf-29). 

Therefore, when the system size is increased beyond this length, the size 
of the excited clusters should stay constant independent of the system size. 
Some results of computation along this line of consideration were presented 
in They demonstrate the fragility of the ground state in a magnetic 
field. For example, the forced spin-flip typically drags as much as 30 percents 
of all spins at H = 0.2 for L = 12. The size of the excited cluster generally 
increases as the system size increases. However, neither a saturation nor a 
proportionality to L'^ could be identifled clearly. 

Therefore, they tried another constraint where signs of an array of bonds 
along a plane are switched. This corresponds to imposing the anti-periodic 
boundary condition to the system. In the RSB picture, the resulting ex- 
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citation again should extend to the whole system, whereas in the droplet 
scenario, it should not for systems larger than 1{H). They measured the 
probability of the excitation wrapping around the system. It was found 
that the wrapping probability increases for H < 0.6 as the system becomes 
larger whereas it decreases for H > 0.7. Based on this result they suggested 
that He is finite and is close to 0.65 at T = 0. While they did not exclude 
the possibility oi He = considering the finite size corrections that may be 
present, their approach seems promising and it is probably worth putting 
more effort in this direction. 

3.6. Sponge-Like Excitations 

While the droplet picture presents a clear "real-space" image for the low- 
lying excitations, we have not discussed how low- lying excitations should 
look in the RSB picture. The authors of discussed geometrical prop- 
erties of low-lying excitations that are compatible with the RSB picture. 
They proposed that low- lying excitations may have a "sponge" -like shape. 
A sponge-like cluster has a characteristic length scale Ic- They argued that 
excitations smaller than this scale obey the droplet scaling whereas be- 
yond this scale the droplet predictions do not apply. A sponge-like object 
is defined by the following properties. 

(1) It is a connected object, and so is its complement. 

(2) It occupies a finite fraction of the whole system's volume. 

(3) It spans the whole system, and so does its complement. 

(4) It has a characteristic length scale different from the system-size and 
the lattice constant. When coarse-grained beyond this length scale, 
the object occupies the space uniformly whereas below this scale non- 
uniformity can be seen. 

(5) Its surface contains a finite fraction of its volume, i.e., the surface-to- 
volume ratio is finite and independent of the system size up to a finite 
size correction. 

They argued that in finite dimensions, three for instance, there may be 
excitations of sponge-like clusters with excitation energy of 0(1), indepen- 
dent of the system size. At the same time, the energy barrier that must 
be overcome in order to excite such a cluster diverges as the system size 
increases. Namely, two pure states in the RSB picture can be transformed 
into one another by flipping spins in one or more of the sponge-like clusters. 
Since they occupy a finite fraction of the whole system, as soon as thermal 
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fluctuations are introduced, they give rise to the non-trivial structure in 
the overlap distribution function, P{q). In addition, because of the prop- 
erty (5), such excitations yield a non-trivial functional form for the link 
overlap distribution P{qi). 

Obviously, the last two properties (4) and (5) above are closely related 
to each other, because the surface-to-volume ratio may be determined by 
the characteristic length scale as S/V ~ Z^^. It is important to estimate 
the length scale, /c, in order to check the applicability of the picture to 
individual cases. An attempt was mad to estimate this ratio numerically 
by considering the EA spin glass model with a Gaussian bond distribution. 
For each bond realization, the ground state was obtained using a heuristic 
optimization technique. Then, the ground state search was performed once 
more, this time under the constraint that a pair of spins (chosen randomly) 
have the opposite relative orientation as compared to the one in the first 
ground state. 

In the droplet picture, the difference of two ground states is supposed 
to be a droplet of a size that is smaller than the distance of the two chosen 
spins, whereas in the sponge picture it should be identified with a sponge 
if one of the chosen spin is included in the sponge and the other is not. In 
the sponge picture, such an event should happen with a pfinite and system- 
size-independent probability because at the scale larger than Ic the sponge 
is uniform and the fraction of the sponge is system-size-independent. 

Calculations of the system up to L = 11 were performed.!^ In order to 
check the property (3), the authors measured the probability of the event 
that the excited cluster spans the whole system. It was found that the 
probability does not show a strong system size dependence, indicating that 
large excitations with 0(1) excitation energy exists. 

3.7. TNT Picture — Introduction of a New Scaling Length 

However, their other findings^^ indicated that the sponge picture defined 
above does not exactly describe the model, either. Particularly disturbing 
was the result of the surface-to- volume ratio of the system-spanning clusters 
obtained by the procedure mentioned in subsection 13. 61 the ratio decreases 
as a function of the system size. A power-law yielded a reasonable fit, i.e. 
S/V oc L"" '^, suggesting Ic (defined as the V/S) tends to infinity in the 
limit L — > oo. If one compares the ground state with the (sponge) excited 
state and focus on a small box located around the center of the system, the 
domain boundary S would never pass through this box in the infinite system 
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size limit. The link-overlap is therefore always complete (i.e., 1) within the 
box. This is in agreement with the result®^ mentioned in subsection 13.41 
The increase of the characteristic length l^, was also confirmed in |1U(J| . It 
follows that the distribution of the link overlap should have a trivial delta- 
peak structure while that of the ordinary overlap may have a non-trivial 
structure. This new scenario with the new system-size-dependent scaling 
length lc{L) is callecP the 'TNT' scenario, an abbreviation of "Trivial 
P{qi) and Non-Trivial P(g)". 

A closely related calculation was done by the authors of |1(J1| , who per- 
formed a Monte Carlo simulation of the EA model with the Gaussian bond 
distribution. They reached rather low temperatures (T ~ 0.2 J) employing 
the exchange Monte Carlo method.EH The overlap distribution P(g) was 
measured. The zero-overlap probability P(0) stays constant, in contrast to 
the prediction of the droplet picture (P(0) cx L^^ with 9 ~ 0.2), consistent 
with the RSB picture and the TNT picture. However, the temperature de- 
pendence of P(0) at a fixed system size agreed with the droplet picture, 
i.e., P(0) cx T. They also obtained the link-overlap distribution. Its width 
turned out to be proportional to where was estimated as 

/i; = 0.76(3) 

from an extrapolation to T = 0. The TNT picture discussed above assumes 
the presence of excitations with size- independent energy, i.e., 9 — Q whereas 
the surface-to-volume ratio of excited clusters goes to zero. Based on this 
picture together with the above estimate of /z/, they obtained (by setting 
61 = in ©) 

d-ds = 0.38(2) 

which is consistent with the zero-temperature calculatiorP^d— dg = 0.42(2). 

3.8. Arguments Supporting the Droplet Picture 

In jl()2| , the notion of pure states was re-examined for the disordered sys- 
tems. The relationship between the appearance of Pj{q) and the num- 
ber of pure states was discussed. It was suggested that Pj{q) is an er- 
roneous indicator of the multiplicity of pure states. In the ferromagnetic 
two-dimensional Ising model below the critical temperature, one has an 
example of P{q) being non-trivial while there are only two pure states. 
Namely, when the anti-periodic boundary condition is imposed in both the 
directions, because of the arbitrariness of the position of domain bound- 
aries, P{q) includes a continuous part. On the other hand, in the three- 
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dimensional random field model, one can see the case where Pj{q) consists 
of a single delta peak although there are two pure states. This is because the 
free-energy difference between the two pure states, one with (5*^) — m > 
and the other with (Si) — — m, diverges being proportional to N^^"^. There- 
fore Pj{q) for a given sample consists of only one delta peak. 

The argument was elaborated in j2Hl()) , in which the authors exploited 
the translation ergodicity. They proved that any translationally invariant 
quantity measured for a particular bond realization J is equal to its bond 
configuration average, provided that the distribution of each bond is in- 
dependent and spatially uniform. Since Pj{q) is obviously translationally 
invariant, it follows that Pj{q) is self-averaging, i.e., does not depend on J. 
This result is in contrast with what we know for the SK model. Further- 
more, they argued that since Pj{q) is self-averaging it would be improbable 
that P.j{q) has a continuous part. The reason for this is that there are only 
countably many pure states and therefore P{q) consists of countably many 
delta peaks. Consequently, we would have to choose many but only count- 
ably many numbers from the interval — 1 < q < 1 to locate these delta 
functions. They argued that existence of such countably many "preferred" 
locations, yet independent of J, are very implausible. Their argument made 
it clear that we must be careful about using the notion of the pure state in 
the disordered systems. 

However, it was pointed ouPSlthat the very existence of the pure states 
that was implicitly assumed in |21) is questionable. This means that the pure 
states that appear in Q may have only a metaphoric meaning, and should 
not be taken too strictly for disordered systems. From this point of view, 
one may say that the non-existence of the unique thermodynamic limit, 
or, more specifically, a chaotic system-size dependence of various quantities 
may signify the RSB nature of the spin glass systems. 

4. Models in Four or Higher Dimensions 

Numerical simulations are easier in four or higher dimensions than in three. 
This is presumably because d = 4 is well separated from the lower critical 
dimension. The existence of a finite-temperature phase transition in four 
dimensions was established already in an early numerical study where a 
clear crossing in the Binder parameter, defined in could be observed. 
The critical point was located at Tc = 1.75(5) and the critical exponents 
were estimated to be z/ = 0.8 ± 0.15 and ry = —0.3 ± 0.15 for the Gaussian 
bond distribution. The scaling relation ^/v — 2 — 7] yields 7 — 1.8(4) 
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which agrees well with the result 7 = 2.0(4) of a high-temperature series 
expansion^°* for the ± J bond distribution. More precise estimates of the 
critical indices were obtained in |1(J5| : 



Based on an off-equilibrium simulation, they also obtained an estimate of 
the EA order parameter below the critical point as a function of the tem- 
perature. A non-vanishing order parameter excludes the possibility that 
the transition is of the BKT-type. As for the nature of the low-temperature 
phase, the observation made in jl()f)j that P(0) is system-size independent 
below the critical temperature indicates that the droplet picture is inap- 
propriate for the four-dimensional spin glass. For a detailed discussion on 
other results obtained before 1997, the readers are referred to the review 
|TO3|. Here we only mention one of the latest computations for clarifying 
the nature of the low-temperature phase in four dimensions. 

The replica exchange method was used in |1U1) to study the 4c? spin 
glass model with a Gaussian bond distribution at rather low temperatures 
T ^ 0.2 J ~ O.lTc. The zero-overlap probability P(0) was measured in 
four dimensions as well as in three. In both cases P{0) was found to be 
independent of the system size for a fixed temperature, and a roughly lin- 
early dependent on the temperature for a fixed system size. The authors 
of |101| summarized their results as being consistent with the TNT pic- 
ture mentioned above. They also estimated the exponent characterizing 
the size dependence of the variance of the link-overlap distribution, namely, 
6 + 2(d~ dg) where dg is the fractal dimension of droplet boundaries (see the 
discussion in subsection 13.4(1 . The extrapolated value to the zero tempera- 
ture was 9 + 2{d — ds) = 0.35(6) while they did not rule out the possibility 
of this value being zero. 

For the ground state properties in four dimensions, we refer the reader 
to |107j , in which the ground states of the model in four dimensions up to 
L = 7 were computed. The stiffness exponent was estimated as 9s = 0.64(5). 

Although somewhat misplaced in this section on higher dimensional 
models we mention here the one-dimensional spin glass model with long- 
range interactions, i.e. Jij = Ca^ij/r^, where ttj is a Gaussian random 
number with zero mean and variance one and the Euclidean distance 
between spin i and j on a finite ring (of perimeter L) embedded in the 
two-dimensional space. This system, on the borderline between finite- 
dimensional and mean field models, was studied in |1()8I1()9| for system 
sizes up to L = 512. Finite temperature simulations as well as the study of 



Te = 1.80(1), 1/ = 0.9(1), 1] 



0.35(5), and 7 = 2.1(2). 
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low-energy excitations yielded results that are inconsistent with the droplet 
picture but (partially) consistent with the TNT picture. 

5. Aging 

The out-of-equilibrium dynamics of spin glasses has become a very rich 
field in the recent years and an excellent review on the intensive theoretical 
work that has been performed on it until 1998 can be found in |110| and 
an overview over the experimental situation until that date can be found in 
|111I112| . Since then a number of interesting developments have occurred 
and we will focus on them with the prerequisites necessary to understand 
them. We start with the theoretical concept of a length scale that evolves 
in time and is flexible enough to account for a number of numerical and ex- 
perimental results. Then we focus on two-time quantities that are typically 
measured experimentally and their behavior in various temperature proto- 
cols (i.e. aging histories) — showing effects like memory and rejuvenation. 
Finally we discuss the theory for violations of the fluctuation-dissipation 
theorem and the numerical and experimental evidences for it. 

5.1. A Growing Length Scale During Aging? 

Spin glasses have an order parameter which is (cum grano salis, i.e. dis- 
regarding the potential complications arising from a possible RSB) qea = 
[{Si}'^]. In comparison to structural glasses^^'^ this is a very lucky situation, 
not only as a starting point for equilibrium theory, which we have discussed 
in last sections, but also for non-equilibrium dynamics. An example is the 
common picture of the dynamical evolution of a system out of equilibrium 
quenched into the phase with non- vanishing qea, in which larger and larger 
regions of space become ordered.^^ Thus it appears natural to postulate 
a time-dependent length-scale L{t) for these equilibrated domains for spin 
glasses. We choose such a picture, motivated by coarsening systems,'^^^ as 
our starting point for reviewing the out-of-equilibrium dynamics of spin 
glasses. 

Within the droplet theorjW^ the slow out-of-equilibrium dynamics of 
a finite-dimensional spin glass is due to thermally activated growth of lo- 
cally equilibrated regions. The growth of such regions is supposed to hap- 
pen through domain wall movements as in the context of pinned domain 
walls in random field systemP^^ or in elastic manifolds in a disordered 
environmenl^^ via a thermally activated process overcoming a free energy 
barrier. It is assumed that the typical energy barrier B^i^t) scales with the 
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typical size L[t) of the domains reached after a time t as 

^ A(L(t)/io)'^, (14) 

where tp is the barrier exponent characteristic of the particular system under 
consideration and A is some constant. Since the dynamics is activated, the 
typical time to overcome a barrier B grows exponentially with B /T^ T being 
the temperature; ^l, activated ~ Gxp{BL/kBT). Hence one would expect that 
after a time t domains of size 

L{t) - {-^Ht/ro)) (15) 

are equilibrated (tq being a microscopic time scale), i.e. the typical domain 
size grows logarithmically slowly with t. Numerically one can determine 
the typical domain size L(t) directly via the spatial two-replica correlation 
function 

1 ^ 

i=l 

where a and b denote two replicas of the same system (i.e. the same dis- 
order realization). For the first time this spatial correlation function was 
studied numerically for the two-dimensional EA spin glass model in |118| . 
and for the 3d EA model in pl9 120 12H, later also in 122 1 2Hll24ll25j . 
The extraction of the typical time-dependent domain-size, L{t), from the 
correlation function G(r, t) given in 1)16(1 is itself a delicate issue. In the first 
studies in |119I118I12U| the definition Lxit) = J dr Gxir, t) was used where 
the subscript T signifies the temperature dependence. This form would be 
justified if Grir^ t) is of a pure exponential form Grir, t) ~ exp(— r/LT(i)). 
Also a scaling form GT{r,t) ~ g{r / LT{t)) has been checked in these early 
works using the length scale obtained via the integral method and a good 
data collapse was obtained for all parameter values (times and temper- 
atures) used. A more flexible functional form was assumed in jl21j that 
was also capable of fitting more accurate estimates of Gt {r, t) obtained 
laterPmni 

GT{r,t)^r-P^^'>g{r/LT{t)). (17) 

The best-fit values of temperature exponent (3(T) were found to be 
constant around 0.5 in three dimensions and decreasing from 1.6 to 0.9 
in the temperature range 1.0 to 0.5 in four dimensions (where Tc « 1.8). 
In |125| only the tail of GT{r,t) was used to extract Lrit) assuming a 
pure exponential form for it. All in all, these different forms do not cause 
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large variations in the estimates of Lxit), the reason simply being that 
the typical domain sizes reached in the times t accessible to Monte-Carlo 
studies to date are only a few lattice spacings. It should be noted that 
these functional forms imply that Grir, t) decays to zero with the distance 
r even in the limit t oo, whereas in the droplet theory one would expect 
a non-vanishing large distance limit: lim^^oo limt^oo GT{r,t) = q%A{T) + 
0{r^^). Functional forms respecting such an asymptotic behavior can be 
also devised to fit the data of Gt {r, t) wellP^ 

Suppose that we have now estimates of L(t) obtained in one way or 
the other from numerical simulations or from experiments. (Experimen- 
tally there is no straightforward way, by which one could possibly measure 
G(r,t) or its space and/or time Fourier transform via scattering or simi- 
lar techniques. Still there is an interesting experimental development^^ on 
this point which we discuss further below.) When fitted to the logarithmic 
growth law (|15|l the estimated value of "0 was approximately 0.1^221 and was 
only acceptable if a finite offset length was introduced. However, a power 
law form with a temperature dependent exponent 

L{t) cx r(^) (18) 

fits the data well, with a{T) « 0.16 T/T^. ■EOUaU This observation may indi- 
cate that we must replace the power-law length-dependence of the barriers 
(I14II by the logarithmic dependence 

cx Alni(i) (19) 

as suggested in |128ll29j . or it may serve as a motivation to modify the 
simple domain growth picture discussed so far in a way that has been first 
discussed in |124| , later also in |1H()I181I12()I125| . The essential idea in these 
discussions is that the initial coarsening process is still influenced by critical 
fluctuations — which is plausible since (i) the temperature that are usually 
studied are not too far from the critical temperature T/Tc > 0.7 and (ii) in 
the initial stage the equilibrated length scales are still very small. 

Close to Tc the equilibrium correlation length is given by ^{T) ^ 
Lo\l — T/Tc\~^, and even at T < Tc the dynamics on length scales L < £^ 
might be dominated by critical fluctuations rather than activated processes. 
Hence one might expect that as long as i < ^ the typical size of equili- 
brated domains grows with time t as L{t) ~ /o(i/io)^^^ (^o and to being 
microscopic length and time scales, respectively) and only for L{t) > ^ 
activated dynamics obeying H15|l sets in. Thus a crossover from critical to 
activated dynamics happens at a time of order To{T)/to ~ (^(T)/Zo)^, when 
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the typical domain size reaches C(T): 

L{t)/i{T)^g{t/T^{T)) with g{x) 



(20) 



^i/z for a; ^ 1 
(Ina;)!/'^ for a; > 1 

Thus, if one plots L{t) ■ [Tc - TY versus t ■ (Tc - TY" one expects a data 
collapse for the time dependent typical domain size at different tempera- 
tures. Indeed such a data collapse has been observed in numerical studies 
of the 4d EA modelpa 
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Fig. 6. Left: Scaling plot of L(i) of the id EA model according to 1201 . where the 
SG transition temperature Tc = 2.0 and the critical exponent v = 0.93 are fixed, but 
the dynamical exponent z is obtained to be 4.98(5) from the best scaling. (From 11251 .) 
Right: Growth laws of the coherence length in the 'id EA model with fitting curves 
according to EJ. (From IT^ .) 



One can go one step furthei^^^^ and write down a more explicit form 
for the relation between length and time scales, similar to \ib\ but now 
taking into account critical fluctuations: 

/A(T)L''^' 



t{L) K, tqL^ exp 



with A(T) w Ao/^{TY (21) 



V kBT 

where Aq is an energy scale of order Tc- Here, following |115| . the tempera- 
ture dependent free energy scale A(T) is assumed to vanish at Tc since one 
might expect that, in analogy to random bond ferromagnets, the surface 
tension of the domain walls vanishes at Tc- The pre-factor tqL^ reflects the 
critical dynamics on short length scales (when the exponential term is still 
small) and might be interpreted as a renormalized microscopic time scale 
To due to critical fluctuations. In this way, we obtain t{L)/fo « f{L/S^{TY, 
with f{x) = exp{cx^) (and c = Ao/fcsT) in reminiscence of (jSOJ. Note, 
however, that now the short time critical dynamics is absorbed into tq. 
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In |126| numerical data for the typical domain size L{t) for the 3d and 4d 
EA model, obtained via Monte-Carlo simulations, were successfully fitted 
to the form For the barrier exponent, ip — 1.0 was obtained in three 
dimensions using ly = 1.65 and z = 7.0, while in four dimensions tJj — 2.3 
was obtained using i/ = 0.8 and z = 5.9. 

As we mentioned above, experimentally it appears to be impossible 
to have a direct access to spatial correlations in spin glass models. How- 
ever an indirect way to estimate a typical size of correlated volumes has 
recently been suggested in |127j . When measuring the waiting-time (i^,) 
dependent thermo-remnant magnetization MTRM(i, ^tu), which is a two- 
time quantity discussed further below, its logarithmic derivative S{t) = 
—dMTKM{t,tw)/d{lnt) at the time t after the field has been switched off 
has a characteristic peak at an effective waiting time t°f . One observes that 
the peak position depends on the strength of the field that is applied. Since 
this time scale is related, via the thermal activation Trdax ~ '''o exp(A/fcBT), 
to the typical value of the free energy barriers that can be explored over 
experimental time scales: A{tw) ~ kBThi{tw/To), the shift in the peak posi- 
tion of S{t) contains information on the field dependence of the free energy 
barriers : A{t^)~E^ kBT\n{t^ /tq), where in [HT) it is assumed that 
is the the magnetic (or Zeeman) energy associated with the change in the 
magnetic field. The latter is connected to the number Ng of spins that are 
involved in the barrier shift induced by the field change via ~ NgXicH'^, 
where Xic is the field cooled susceptibility per spin and H the magnetic field 
strength. If one assumes that the Ng spins to be effectively locked together 
one gets an estimate for a correlated volume Ng oc i.e an estimate for 
the coherence length ^ as a function of the waiting time 

By analyzing their data for the peak position of S{t) to obtain ^(tu,) 
in the way just described the authors of |127| find that it fits to (|18|l with 
a(T) « 0.169 T/Tg, in agreement with the numerical results reported in 
|120I121| . A fit to the logarithmic form (|15|1 works equally well, but best-fit 
value of ip turns out to be rather large 5), which does not agree with the 
numerical estimates cited above. These experiments probe a time window 
that is different from the one accessible to numerical simulations; each is 
roughly 6 decade wide but centered around times roughly 10 decades apart. 
Thus only the algebraic growth law (|18|l . or the form H21|l that interpolates 
between critical and activated dynamics, would consistently match the data 
of both time windows. 
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5.2. Two Time Quantities: Isothermal Aging 

The out-of-equilibrium properties of glassy materials, quenched rapidly be- 
low the glass transition temperature Tg and then aged isothermally (i.e. at 
constant temperature T < Tg), manifest themselves most prominently in 
two-time quantities, typically response functions, susceptibilities or corre- 
lation functions. The magnetic response function is defined as 



(22) 

hi=0 



In experiments as well as numerical simulations one usually does not apply 
field pulses after a waiting time but one switches a constant magnetic field 
on or off after a waiting time — the corresponding susceptibility is then 
related to the response function via a time integral 

Xit + tw,ty,) ^ J dt' R{t + t^,t') (23) 

If one switches the field on at time after the quench and measure the 
susceptibility, one calls it zero- field- cooled (ZFC), whereas if it is switched 
off it is called field-cooled (FC). In the former case the magnetization, which 
is within the linear response regime simply related to the susceptibility via 
M[t, tw) — h-x{t, iu)), increases with the time t spent in the field, whereas in 
the latter case it decreases with the time t after the field has been switched 
off — it is then also called the thermo-remnant magnetization (TRM). 

Very useful for the development of theoretical concepts is the two-time 
spin autocorrelation function 

C(t + t^,t^) = A^-i^(5,;(i + t^)^,(i„)) . (24) 

i 

In equilibrium it is related to the response function, the susceptibility and 
the ZFC/FC magnetization via a fluctuation-dissipation theorem (FDT). 
(We discuss the violation of the FDT in the subsection l5.4n However, since 
we are out of equilibrium, we have first to regard C as being an independent 
quantity. 

The huge amount of experimental data that have been collected are 
nicely over- viewed in |llllll2j . Let us here simply state the main results. 
For the field cooled magnetization one observes that the remnant magneti- 
zation can be decomposed into a stationary part Mcq{t) that depends on t 
only and an aging part Maging(i -|- i^,). Sometimes this decomposition is 
assumed to be multiplicative, M(t-\-tuj, tw) = -Meq (i)' -Paging (i+^u;, ^u;), as in 
the universal short time dynamics for coarsening at the critical point,^^^^^ 
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and sometimes additive, M{t + t^,tw) = M^qit) + Maging(i + t^.t^), as in 
coarsening dynamics of pure system^^ and in the out-of-equihbrium the- 
ory of mean- field spin glasses Since Maging(i + tw^ttu) is approximately 
constant for t <^ both forms usually fit the data well, although they 
are fundamentally different. For instance the multiplicative decomposition 
is incompatible with the existence of a plateau in C{t + tii,,tw) that occurs 
in mean-field theories of spin glasses at low temperatures for long wait- 
ing timesr'""'^ Analogously it is also hardly compatible with a non- vanishing 
long time limit of limt^^co C(t+tw, tw), which is, within the droplet theory, 
identical to the order-parameter qea- 

The most remarkable feature of the aging part is that it scales to a good 
approximation with t/t^ 

-Paging ~^ tw T tw) ~ Afaging (^/^ii> ) (25) 

In ZFC experiments, instead of applying a constant field after a waiting time 
tu,, one usually applies a small oscillating field with frequency uj, which 
essentially means that one measures the t- Fourier transform xi'-^yt^) of 
x{t + tw,tw)- The frequency uj plays the role of an inverse observation time 
1/t, with t <^ t^. The waiting time is usually in the range of hours 
or more, while cj is in the range of O.lHz to lOOHz, i.e., t is less than ten 
seconds. Again x{t~^t^,t^) can be decomposed into a stationary part Xcq{<^) 
and an aging part Xaging('-^, ^to) that now scales to a good approximation 
with Lotw'- 

Xaging ("^ : i tu ) ~ Xaging(^ito) (26) 

corresponding to the aforementioned t/tw scaling for Mpc- Numerically 
one also has access to the spin autocorrelation function H24|) . which can 
also be decomposed into two parts. In |120| it has been shown that (i) the 
equilibrium part Ccq(i) decays algebraically with a very small temperature- 
dependent exponent x{T) and (ii) that Cago(i+ito, ^to) again obeys to a good 
approximation the t/t^ scaling: 

Ccq(0 « ^i"""^^^ and Cagc(^ + i»,t«,) ~ S(iA») , (27) 

where the scaling function c behaves as c{x) (x x~^'^'^^ with A(r) being a 
temperature-dependent exponent much larger than x{T). This behavior was 
also found for the 4d EA model. Later studie^^^ focused on small but 
systematic deviations in the aging part from a simple t/t^ scaling behavior; 
the values of Cage {t + tw^tw) for fixed ratio t/t^ show a slight tendency to 
decrease with increasing t^^ which is called sub-aging and is interpreted in 
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terms of an effective relaxation time t^ff that is smaller than the actual 
waiting time tw We shall return to this point later. 




10"^ 10"'* 10"^ 10"^ 10"^ 10° 10^ 10^ 



t/tw 

Fig. 7. Left: Aging part of the FC magnetization Mpcit+twtt-w) measured in an AgMn 
spin glass: the estimated stationary contribution Meq(t) has been subtracted from the full 
measured value. The data is plotted versus t/t^. The deviations from perfect t/tm scaling 
occur for t ^ tw (sub-aging) and can be accounted for by choosing h{t + tw) / h{tw) with 
h{t) given in 1301 as a scaling variable. (From lllll .') Right: Autocorrelation function 
C{t + tw,t^), devided by t"^*^) oc Coq(t), versus t/tw for T = 0.4 in the 3d EA model. 
(From inni-) 

Within a domain growth (or droplet) picture of the out-of-equilibrium 
dynamics of spin glasses only one characteristic length scale L{t) for a 
given time t after the quench is assumed to exist. One would then expect 
that two-time quantities depend only on the ratio of the two length scales 
that are equilibrated within the two times t„, and t + t^: 

0,sing{t + t^,t^) ^d{L{t + t^)/L{t^)) , (28) 

where 0(t + tw, tw) is any two-time observable (such as R, X: M or C) and o 
the corresponding scaling function. The first observation is that an algebraic 
growth law, 1)18(1 for instance, would then be completely compatible with 
the t/tw scaling of the aging part of various observable reported so far. Note 
also that as long as L{t) depends algebraically on t the use of the scaling 
variable L{tw) / L{t) is equivalent to using L{t + tw) / L{tw), both end up to 
be some function of t/ tw- 
in the case of logarithmic domain growth on the other hand, 
the scaling with L{tw)/L{t) is different from the scaling with L{t + 
tw)/L{tw), and one should note that droplet scaling,^^ as it is usually 
applied to the experimental and numerical data, is assumed to work with 
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L(t^)/L(i)lSniSlSIl where L(t^) is the typical size of the domains after 
waiting time and L{t) is the typical size of the droplets being activated 
in the time t or being polarized by an ac field of frequency lu = 1/t. The 
basic physical idea behind this scaling form is that the presence of a domain 
wall effectively reduces the excitation gap of droplets that it touches. 

However, assuming the simple logarithmic growth law H15|l . L(t) « 
((/cbT/A) ln(t/To))^/''', the numerical as well as the experimental data turn 
out not to scale with L{t^)/L{t) as suggested by the droplet theor j^SGSl (it 
does not scale with L(t + 1^)/ L{tw), either, which is incompatible with the 
experimental observation of sub-aging) . For some time this observation was 
taken as an indication that the droplet theory might be inappropriate to de- 
scribe the out-of-equilibrium properties of spin glasses, at least for the time 
scales accessible to experiments and to numerical studies. However, recently 
it was pointed out^'^^ that, taking into account the crossover from critical 
to activated dynamics already discussed in the last subsection, the appar- 
ent inconsistency of the data with the droplet theory disappears. There 
experimental data for the (imaginary part of the) susceptibility x(w,iu,), 
measured for the 3d Ising spin glass Feo.sMno.sTiOa and the 3d Heisenberg 
spin glass Ag(ll% at Mn), were presented. The aging part Xaging(^' ^m-) 
(obtained after subtracting a fitted equilibrium part Xeq) was shown to 
scale as L{tw) / L{1 / lo) with L{t) following the logarithmic growth law H15|l. 
where the microscopic time and energy scale tq and A replaced by a tem- 
perature dependent characteristic time and energy scale, Tc(T) and A(r). 
This replacement captures the features of the critical, short-time dynamics 
close to the critical point: Tc{T) - r^^^ and A(r) - ^"'^ (^ = \l-T/Tc\^'' 
being the correlation length) , which is in the same spirit as (I2n|l and H21(l . 

The authors of |131| found a very good data collapse when plot- 
ting their data for (x"(w,tM,) — Xcq)/x"('J-', i^) versus the scaling variable 
ln(t/Tc(T))/ln(w-irc(T)) (which is L{t) / L{1 / uj) , since A(r) and ksT can- 
cel in this ratio). Given the critical exponents v and z from earlier studies, 
the fitting parameters are ij:, the whole equilibrium part Xeq (for which an 
explicit functional form predicted by the droplet theory was used) and a mi- 
croscopic time scale Tm- They found "0 «:! 1.9 for the Ising spin glass, which 
was much larger than what had been reported before in experiments^^® (see 
also [13011321 1 and in the numerical simulations discussed in the previous 
subsection, and tp ~ 1.3 for the Heisenberg spin glass. The authors of |131| 
showed that their experimental data are also compatible to the Lut„ scaling. 

So far, we have discussed the compatibility of numerical and experimen- 
tal data with a simple (or simplified) domain growth picture. Mean field 
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models, however, show a more comphcated behavior fSMSI The mean field 
theory for the out-of-equilibrium dynamics predicts "ultra-metric behavior" 
in the time domainPHSnEM 

Paging H" t 

i 

where the infinite sum over the index i refers to various large time 
sectorP^ESl defined by the ratio of hi{t + tw)/hi{tw) being of order 1, 
where the different (unknown) functions hi{t) represent different physical 
mechanisms at work during aging and need not necessarily to be related 
to domain sizes. They are monotonously increasing and grow differently 
such that when < hi{t + tw)/hi{tyj) < 1 holds for an index i then 
hj{t + tw)/hj{tw) = 1 for all larger indices j > i — implying that if two 
times ti and t2 belong to the time sector defined by hi and ^2 and ^3 to 
the one defined by hj with j > i, then ti and ^3 also belong to the sector 
defined by hi. One simple example for hi{t) is 

h,it) = exp{(tAo)'-^7(i - ^^^)} , (30) 

with < fii < 1 {fi ~ 1 yields a i/i^-scaling, /i = time-translational 
invariance. 

A consequence of such an Ansatz is for instance a hierarchy in these 
large-time sectors, called dynamic ultrametricity, which for the correlation 
function C means that in the limit of large times ti < t2 < 

Cih,h) = min{C(t2, h), C{h, h)} (31) 

holds. A number of consequences can be drawn from this Ansatz.^^^^ 
Among them is a particular form of the violation of the fluctuation- 
dissipation relation, which we discuss later. Let us here just state that mean- 
field theory predicts a richer scenario for glassy out-of-equilibrium dynam- 
ics than the domain growth picture with only one waiting-time-dependent 
length-scale. However, as rich as this scenario is, all experimental data for 
real spin glasses and all numerical data for spin glasses with short-range 
interactions obtained for isothermal aging up to now can, to our knowledge, 
be scaled nicely with one single large-time domain, i.e. with only one term 
appearing in the infinite sum H29(l . 

For instance, the aforementioned sub-aging property of two-time quan- 
tities that has been observed in experimental datc^^^^ (and also in an- 
alytically tractable coarsening, non-spin glass, modelP^ES) can easily be 
accounted for by taking only a single term in the sum 

with h{t) being 

the function defined in and with the exponent /z now being a fitting 
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parameter. This form was actually used in the earliest experiments on aging 
in polymer glasses/**^ and later also in spin glass experiments^^ The expo- 
nent fi commonly turns out to be not very different from 1, accounting for 
the fact that the deviations from simple t/tw scaling are usually very small. 
In numerical simulations of 3d models the scaling with t/t^, i.e. fi = 1, 
is almost perfect. In four dimensions the numerical data for C{t + t^^t^) 
tend to show super-aging^^ i.e. the values of Cage + iu, , tu, ) for fixed ra- 
tio t/t^} show a slight tendency to increase with increasing t^. This can 
be interpreted by an effective relaxation time that grows faster than i^,, 
and it can be shown that the aging part Cagc(i + tw,tw) scales nicely with 
h{t + tyj) / h{tw) , when /i > 1 for hit) in from (|30|l is chosen (which bears 
unfortunately mathematical inconsistencies^). However, the data can also 
be scaled with a domain growth scaling form H21(l and a non- vanishing long 
time limit for the equilibrium part.^^ 

5.3. More Complicated Temperature Protocols 

From what we learned from isothermal-aging experiments described in the 
last subsection, it follows that the dynamics of a spin glass depends crucially 
on the age of the system, i.e. the time t^ that it spent in the glass phase, i.e. 
at T < Tg. Moreover, as experiments have impressively demonstrated,^Sili2I 
it depends sensitively on changes of temperature (and also field) during this 
waiting time tw A major goal of their systematic study described in this 
subsection is to understand the different phenomena that are measurable 
by applying different protocols (i.e. temperature variations during aging) 
and relating them to a theoretical picture, if not a theory, of the out-of- 
equilibrium dynamics of spin glasses. 

The simplest experiments in this direction are small-temperature-shift 
experiments.^^ The protocol is as follows. First the system is rapidly 
quenched below Tc to a temperature Ti where it is aged for a waiting time 
tw Then the temperature is shifted to a new temperature T2 = TiztAT and 
the measurement is started. The first systematic numerical study of such a 
protocol was performed for the Sd EA model in |147j and a little later in 
[126j . It was observed that the decay of the correlation function is slower 
for negative shifts (T2 < Ti) when compared with the same function aged 
isothermally at Ti. For small shifts AT the functional form of C{t + tw,tw) 
is the same as that for isothermal aging, and it can be matched with it 
using an effective waiting time < tw for T2 < Ti. If one assumes that 
the barriers that the system can surmount during a time tw at temperature 
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Ti are of the same size as those that it can surmount during a time <^ 
at temperature T2, one obtains a good agreement with the numerical data, 
only if one replaces the microscopic time tq by the typical time scale for 
critical dynamics as in 1)2 1(1 

Obviously such an interpretation of the numerical results implies that 
successive aging at two different temperatures will add to each other in a 
full accumulative way. The effective age is a monotonically increasing 
function oitw and also depends on the pair (Ti,T2): ij^f = fTi,T2{tw), and 
tyj can be related to by the inverse: t^^ = /j^^^y^ (i^f )■ In |148| twin- 
experiments — negative shift from Ti to T2 and positive shift from T2 to 
Ti — were considered. A criterion for purely accumulative aging would be 
^TiT2 ~ /t^.Ti {t) and in contrast to the aforementioned numerical results 
deviations from it were observed in these experiments for AT/Tc > 0.01, 
i.e. still very small shifts. The authors then concluded that positive and 
negative temperature shifts both cause a restart of aging and these findings 
were then interpreted as a symmetrical temperature-chaos effect .^^1^^^^ 

The temperature chaos in spin glasses is one of the basic ingredients 
of the droplet theorj^^ and implies that the equilibrium configurations at 
two different temperatures Ti and T2 = Ti± AT below Tg are uncorrelated 
beyond a length scale called the overlap-length, Iat ~ (AT)~^/'», where ^ 
is the chaos exponent. A number of experimental results, in particular ob- 
tained from temperature-shift or temperature-cycling experiments,'^^ have 
been interpreted within this scenario. Numerical estimates of the value of 
C are based on small variations of the couplings Jij (via Jij — > Jij -\- Sij 
with Sij small and random) at T = (^lESI rather than variations of the 
temperature, simply because ground states of spin glasses, in particular 
in d = 2 can be obtained more easily than equilibrated configurations at 
small T. The reported estimates are C « 1 both in two and three dimen- 
sions. While the same exponent was reported for temperature chaos (via 
Monte-Carlo simulations) in the two dimensional case, recent large scale 
numerical studies of the 3d EA model did not show any evidence for tem- 
perature chaos in spin glasses on the length and time scales that could 
be probed.^^ Even if temperature-chaos exists in three dimensional spin 
glasses, the overlap length might be much larger than the length-scales that 
have been equilibrated during a particular waiting time, which then makes 
the chaotic rearrangements due to the temperature shift invisible. However, 
there still might be some effect due to dangerously irrelevant droplets whose 
free energy gaps are quite small, as discussed in 125 150 154 . 

The next temperature protocols that we discuss here are larger shifts 
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and cycles Ti ^ T2 ^ Ti .^Sl xhe first shift experiments was reported in 
|155j whereas the cycle experiments can be found in (1561 . The first quali- 
tative numerical study of this kind was reported in |157| . In the theory of 
the asymptotic out-of-equilibrium dynamics of mean-field spin glass models 
one can understand these effects on the basis of infinitely many time scales 
organized in a hierarchical — ultra-metric — way (which we will discuss in 
section IF^ . 

A systematic numerical study of large-shift and cycle experiments was 
recently performed for the 3d and Ad EA model in |126| . The basic mes- 
sage of large temperature shift experiments is that, independent of the sign 
of Ti — T2, aging is "restarted" at the new temperature. This rejuvena- 
tion effect can nicely be observed in experiments as well as in simulations 
through the measurement of the susceptibility Rejuvenation after 

a negative temperature shift comes from fast modePSEnnS! 

which were 

equilibrated at Ti, but fall out of equilibrium and are slow at T2. Therefore, 
one should expect to see this phenomenon if the spin configurations in the 
equilibrated regions (on length scales < Lxiitw)) are sufficiently different 
at the two temperatures. This mechanism is obviously qualitatively differ- 
ent from the interpretation involving the notion of temperature chaos (see 
above), which implies that length scale smaller than the overlap length are 
essentially unaffected by the temperature shift, while larger length scales 
are completely reshuffled by the shift. In this picture, rejuvenation is thus 
attributed to large length scales and strong rejuvenation effects therefore 
require a very small overlap length. 

No clear rejuvenation effects have ever been observed in simulations of 
the 3d EA niodelP™nTIMl 

This was first attributed to the fact that 
the overlap length was perhaps numerically large, so that no large scale 
reorganization could be observed on the time scale of the simulation. How- 
ever, in the fast-mode mechanism mentioned above, the crucial ingredient 
is the small-scale reorganization due to a temperature shift. When compar- 
ing the spatial correlation function G't(7', itu) at two different temperatures 
and waiting times such that the length scales iTi(iu)i) — iT2(it02)i there is 
nearly no difference even for quite large AT — T1 — T2 in three dimensions 
but significant differences in four dimensions, which is compatible with 
the observation that the exponent (3{T) in H17|l remains roughly constant 
for a large temperature interval in three dimensions whereas it varies signif- 
icantly in four dimensions. This observation suggests that the Ad EA model 
should be more favorable in studying temperature shift/cycle protocols. 
The numerical simulations of the Ad EA model show indeed a clear reju- 
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venation effect in increasing smootfily with AT, wfiich one would 

not expect within the temperature chaos picture. (See FiglHl (Right).) 
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Fig. 8. Left: Series of 'dip' imprinted on the a.c. susceptibility by successive stops 
at different temperatures while the system is cooled. Further cooling 'rejuvenates' the 
system (i.e the susceptibility goes up). However, the dips are one by one remembered 
by the system when heated back. For more details, see 11631 . (The figure is adopted 
from 11321 '). Right: Evolution of the 'a.c' correlation function taken from Monte-Carlo 
simulations of the id EA model in the procedure T = oo — » Ti = 0.9 — ► T2 = 0.4 — > 
Ti — > T2, showing, as in experiments the coexistence of rejuvenation and memory effects 
(from Tm). 



The experimental procedure for cycles is T = 00 — > Ti ^ T2 < Ti — > Ti . 
The time spent at Ti is ts and the time spent at T2 is t'^. The spectacu- 
lar memory effect arises when the temperature is shifted back to Ti. It is 
observed that although aging was fully restarted at T2, the system has a 
strong memory of the previous aging at Ti. The dynamics at Ti proceeds 
almost as if no cycle to T2 has been performed The coexistence of reju- 
venation and memory was made even more impressive in a multiple step, 
or dip experiment in which the temperature was decreased to Ti, then 
maintained for some time, then further decreased to T2 < Ti, again main- 
tained, then further decreased to T3 < T2 etc. At some lowest temperature 
the sample was then re-heated with the same rate as it was cooled before, 
but now without the interrupts at Ti, T2, T3, • • • . As shown in FiglHl (Left), 
on the way down to the final temperature the susceptibility x(a;) as a func- 
tion of the temperature showed clear dips at Ti, T2, Ts, . . ., which is the 
rejuvenation effect. Most remarkably the same dips reappear in xi^) 
the way up, which is the memory effect. The spin glass has memorized the 
individual temperatures at which the constant cooling was interrupted for 
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some time. 

According to |182j the memory effect is a simple consequence of the sep- 
aration of time and length scales, also observable in simpler slowly coars- 
ening systems with activated dynamics.^^ When the system is at T2 < 7i , 
rejuvenation involves very small length scales as compared to the length 
scales involved in the aging at Ti. Thus when the temperature is shifted 
back to Ti, the correlation of length scale LT2{t's) nearly instantaneously 
re-equilibrate at Ti. The memory is just stored in the intermediate length 
scales, between ^^^(i^) and LTiitg). How such a storage mechanism actu- 
ally works in a microscopic model for a finite-dimensional spin glass, like 
the EA model, is a challenging question. Progress has been made for sim- 
pler microscopic models, such as the 2d Mattis model, a simple coarsening 
model,lSl or the directed polymer in a random medium.ESSl Also experi- 
ments on random ferromagnetic systems] ^^^ ^ ^^^ ^ ^^^ ^ ^^^ ^ ^'^^ on random ferro- 
electric system^^^^ and on frustrated systems without a spin glass tran- 
sition at finite temperatureS^^show strong rejuvenation effects and a weak 
memory effect. This indicates that the existence of "many pure states" as 
predicted by Parisi's equilibrium solution of the Sherrington-Kirkpatrick 
model and the interpretation of these non-equilibrium effects based on it 
(in terms of a diffusion in a hierarchical are not necessary pre- 

requisites to observe rejuvenation and memory in the out-of-equilibrium 
dynamics. 

Nevertheless the hierarchical picture might be relevant in more complex 
situations and a concrete implementation was proposed in |175) in terms of 
a thermally activated random energy model (for an overview see Here 
to each level of a hierarchical tree a transition temperature is associated, 
such that for each level the dynamics is stationary for higher temperatures 
and aging for lower temperatures. A small decrease of temperature induces 
some rejuvenation by driving out of equilibrium a new level of the tree, while 
freezing out the dynamics at the upper levels, thereby allowing the memory 
to be conserved. This model was recently studied further,t2SI where it was 
shown numerically that the rejuvenation/memory effect is indeed already 
reproduced with two levels. A real space interpretation of this hierarchical 
tree was proposed in |175ll77j . and further developed in ( 159) in terms of 
a multi-scale dynamics. Low levels of the tree correspond to short wave- 
length modes, which are only frozen at low temperature, while large wave- 
length modes are frozen at a higher temperature and constitute a backbone^ 
where memory is imprinted. In the context of domain walls this picture is 
particularly cleaJ^^ and should apply directly to disordered ferromagnets 



February 2, 2008 5:33 



WSPC/Trim Size: 9in x 6in for Review Volume 



main 



48 N. Kawashima and H. Rieger 

where the slow dynamics comes from the motion of these pinned domain 

^alls ;.166 167 168 169 170 

5.4. Violation of the Fluctuation-Dissipation Theorem 

In a system in equiUbrium the response R to an external magnetic field 
(Eq.(|22J)) and the autocorrelation function C fEa. H24l) ') depend only on the 
time difference t, i.e., R{t + tuntw) = Rcq{t) and C{t + tw,tw) — Coq(i), 
and are related to each other through the fluctuation-dissipation theorem 
(FDT) 

^=.,(0^4^. (32) 

In a system that is out-of-equilibrium this relation is generally not valid 
and the violation of FDT can be parameterized through a violation factor 

X{t + tyj,t^): 

s X{t + tw,tw) dC{t + tw,tw) , s 

R(t + tw,tyj)^ — — , [66) 

where the differentiation is done with respect to , not the time difference 
t as in (|32|l . In analytic studies of mean- field systems it can be showrPSl 
that for large times (t^ oo) X depends on t and only through the 
value of the correlation function. Hence X{t + t^^,t^^) = X[C(i-f t^,, t^,)]. In 
particular, when C > qea it is X = 1 and the FDT is recovered. In essence, 
the asymptotic dependence of X on the values of C alone underlies the 
decomposition into asymptotic time domains discussed in subsection 15.21 
under H29|l . It turns oui^^ that the effective temperature Tcff{t -f tw,tw) = 
T/X{t -\- tw^tyj) is precisely the temperature which would be read on a 
thermometer with response time t (or frequency w ~ 1 /t) when connected 
to the magnetization at time tm. A "fast" thermometer of response time 
t <^ tyj will then probe the stationary regime for which X — \ and thus 
measure the heat-bath temperature. 

Different scenarios of the FDT violation can occur. (A) X{C) ~ 1 for 
all C implies the system is in equilibrium, e.g. in the paramagnetic phase 
of a spin glass. (B) Coarsening systems with qea = m"^ have X{C) = 1 
for C > qea and X(C) = for C < qea- (C) So called discontinuous 
spin glasses, such as the spherical mean field spin glass model with p-spin 
interactions (p > 3),C2Hlhave X{C) = 1 for C > qea and X{C) = xi < 1 
for C < Qea, where xi is a constant. This corresponds to the existence of 
only one asymptotic time domain (see ()29|lV and the two-time quantities 
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are expected to scale with only one (possibly unknown) function h{t^^^ 
(D) So called continuous spin glasses, like the SK model, have X{C) = 1 
for C > qea and X{C') continuously varying, non-constant function of C 
for C<qEA^ 

A convenient way to extract the function X{C) is to measure the auto- 
correlation function C{t + tw^tyj) and the susceptibility x(t -I- i^,, t^,), H23|l . 
for a situation in which a (infinitesimal) small homogeneous external field 
h is switched on after a time t^^, (in which case the induced magnetization 
is zero- field-cooled and xit~^'^w,'tw) — ^v[ah^aMzFc{t + 'tw,tw)/h). In this 
case we have 

Xit + t^,tyj) = / dt'R{t + t^,t') 



= 1 j'^'^ dt' X[C[t + t^,t')]^^^^ . (34) 

As long as the function C{t + t^, t') is monotonously increasing with t' for 
fixed time t + t^, the substitution of X{t + tw, t') by X[C{t + tyj, t')] is legit- 
imate. One should keep in mind that then X[C] still has a tu,-dependence 
(that we could indicate by Xt^[C]), which however vanishes in the limit 
— > CO. Hence one gets (with C{t + i^,, t + 1^,) ~ 1): 

X{t + U,,t^)^ f dCX[C] or for fixed t^: ^ = -^P 

(35) 

which implies that for fixed waiting time the slope of a parametric plot of 
x{t -t- ito, t«,) versus C{t + t^jt^j^) yields Xt^ [C]. In the limit t„, ^ cx) one 
should obtain the desired FDT-violation X[C] ~ limf^^oo -^t„ (C"). 

The first study of the FDT violation function X{C) in the 3d EA model 
was performed in |18U| (the first report of an FDT violation was actually 
already presented earlier in |181j ). The non-constant part of X{C) still 
showed a small but systematic ti„-dependence for the small-C branch of 
the curve X{C). Since no tendency to approach a straight line (scenario 
B or C) could be observed the data were interpreted as an indication for 
scenario D, which implies that the 3d EA model has an out-of-equilibrium 
dynamics that is similar to the SK model. Later studies confirmed this 
result, in three dimensions'^! a,s well as in four dimensions However, 
in jl25| it was proposed that the numerical data for the FDT-violation in 
the 3d EA-model (as well as those for the spatial correlations and the two- 
time quantities) can be interpreted in view of an extended droplet theory 
and that for the 3d EA model scenario B (essentially slow coarsening in 
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C{t,tw) Correlation C(t',t) 



Fig. 9. Left: Parametric plot {T / J)xzFc{t + tni,t^) vs C{t + t^,t^) at T/J = 1.2 
for the Ad EA model. The straight tangent lines represents the FDT. The convergence 
points of the break points of time translational invariance (qeAi {T/ J)xea) a^nd the 
break points of FDT (qo , (T / J)xd) ■ (from EH) Right: FDT-plot obtained from the 
experiments performed in 11851 . Relaxation measurements are plotted versus correlation 
functions for each t' . The dot-dashed line (FDT line) is calculated for T = O.STg = 
13.3-R', from a calibration obtained with a copper sample. The dashed line represents the 
scaling extrapolation for t' — > oo. The branching point with the FDT line, corresponds 
to C = Qea (square symbol, with size giving the error range). In Inset, the same data 
in the whole range, (from flSSl). 

the spirit of the droplet model) is appropriate. (See Fig|51 (Left).) Since 
numerical simulations can cover only a few decades of the waiting time 
tw and there is, as mentioned above, still a significant dependence in 
X{C, tw) in the numerical data,^*'^ ^^'^ a definite statement is hard to make. 

For a clear experimental demonstration of the FDT violation one has 
to measure both; response and autocorrelation function (in jl8f)j possible 
indications of the FDT violation based only on the response have been 
discussed). Such an experimental evidence of the violation of the FDT in 
spin glasses has been presented for the first time in |185j (for a structural 
glass, glycerol, it has been measured in fl87j . for a colloidal glass, laponite, 
in |188j V where autocorrelation function C{t+tjj^, tw) has been measured via 
the magnetization fluctuations in a spin glass sample and the susceptibility 
xit + ti^,tw) in the usual way by applying a magnetic field. (See FiglHI 
(Right).) The x versus C plots show still a strong waiting time dependence 
and give some room for extrapolations to the long time limit — one that is 
suggested by the authors has similarities to the one for the SK-model. On 
the other hand, the aging part of their data for C and x scale perfectly with 
h{tw)/h{t) where h{u) = exp(M^~^/(l — fj,)) with fj, = 0.85, which would 
not correspond to what one expects for the SK model. 
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The action for the long-time asymptotic behavior of the generating func- 
tional for the correlation and response functions is invariant under a time 
re-parameterization t — > h{t), simply because in the long-time regime a 
re-parameterization of time with a monotonous function does not change 
the ultra-metric relation in the time ordering. The global invariance under 
the re-parameterization group (RpG) for mean-field models was recently 
extendecP^^Eini to its local variant for short range models to treat fluc- 
tuations. In analogy to Heisenberg magnets, where the global spin rota- 
tion invariance lead to the existence of spin waves in finite dimensional 
models via the Goldstone modes, the authors proposed that local time re- 
parameterizations play a similar role for the asymptotic dynamics of short- 
range spin glasses. One consequence of this observation is that one expects 
fluctuations in local correlation functions and local susceptibilities 

/h, 

(36) 

where Sr{t) means a slightly coarse-grained spin value at position r and 
time t, i.e. it involves a spatial average of the spin values over a small 
volume centered around r and a time average over a small time window 
centered around t. If the volume is extended to the system size, the usual 
global quantities discussed above are recovered. These spatially fluctuating 
quantities then have a joint probability distribution p{Cr, Xr) that stretches 
along the FDT-relation xCC*) for the global quantities x and C, as has been 
shown numerically in |190| . 

5.5. Hysteresis in Spin Glasses 

Due to the complex energy landscape of spin glasses one expects 
strong hysteresis effects, which were studied recently via Monte-Carlo 
simulation j-'^^^ l ^^^ l ^'^^ l ^^'^ l also using new techniques such as the recently in- 
troduced First Order Reversal Curve (FORC) method!^^^' 

In the zero temperature dynamics of the 2d EA spin glass model 

was simulated with a varying external field H . The procedure used there 
was as follows. The magnetic field is changed in small steps, first downward 
from positive saturation and then upward from a reversal field Hn. After 
each field step, the effective local field hi of each spin Si is calculated: 
hi — 'Y^- JijSj — H. A spin is unstable if hiSi < 0. A randomly chosen 
unstable spin is flipped and then the local fields at neighboring sites are 
updated. This procedure is repeated until all spins become stable. 

The first important observation was of a memory effect in the hysteresis 
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of the 2d EA spin glass model that emerged when the magnetic field was 
first decreased from its saturation value and then increased again from 
some reversal field Hn. It was found that the EA spin glass exhibited a 
singularity at the negative of the reversal field, —Hn, in the form of a kink 
in the magnetization of the reversal curve (See Fig ll()|l . 



1 
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^ 
-0.5 
-1 
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Fig. 10. Reversal curve (solid line) and major hysteresis loop (dotted line) for a two- 
dimensional (2D) EA spin glass with 10^ spins and Hn = —2.28. In the inset a kink is 
seen around -Hr. (From 

One can describe this effect within a phenomenological approach to 
hysteretic systems, the Preisach model.E^In the Preisach model a magnetic 
system is described as a collection of independent two-state (±1) switching 
units, or "hysterons". Unlike Ising spins, which always align with their 
local field, the hysteron's state changes from —1 to -1-1 at a field Hb + He, 
different from the field Hi, — He, required to switch the hysteron from +1 
to —1. Different systems are distinguished by their different distributions 
p{Hb, He) of hysterons of a given bias Hf, and coercivity He- Here p{Hi,, He) 
is the so-called "Preisach function" . 

This function was extracted by a tool developed for analyzing exper- 
imental data of hysteretic systems.'iSSl A family of First Order Reversal 
Curves (FORCs) with different Hr was generated, with M{H, Hn) de- 
noting the resulting magnetization as a function of the applied and re- 
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versal fields. Computing tlie mixed second order derivative p(H, Hji) = 
-(l/2)[d'^M/dHdHn] and changing variables to = {H - Hr)/2 and 
Hi, = {H + H]i)/2, the local coercivity and bias, respectively, yielded the 
"FORC distribution" p{Hi,,Hc)- For phenomenological Preisach models, 
the FORC distribution is equal to the Preisach function. However, the def- 
inition of the FORC distributions is more general, because it is extracted 
from numerical or experimental data, and thus is model- independent. 

Fig. II II shows the FORC diagram of the EA spin glass. The ridge along 
the He axis in the range 1.5 < He < 4.0 corresponds to the kinks of Fig.llUI 
Thus FORC diagrams capture the reversal-field memory effect in the form 
of a ridge along the He axis. In Fig. also the experimentally determined 
FORC diagram of thin films of well-dispersed single-domain magnetic Co- 
7-Fe2 03 particles provided by Kodak Inc is shown. It clearly exhibits the 
horizontal ridge associated with the reversal- field memory effect. This strik- 
ing similarity between the experimentally determined FORC diagram of 
the Co-7-Fe203 films and the numerically determined FORC diagram of 
the EASG indicates not only that Co-7-Fe203 films exhibit reversal-field 
memory but also that frustration may be a component of the physics of the 
Co-7-Fe203 films. Note that for instance the 2d random field Ising model 
does NOT exhibit a reversal-field memory and that its FORC diagram has 
a vertical ridge rather than a horizontal one.^^ 




Fig. 11. Left: FORC Diagram of the EASG. Note the ridge along the He axis. Right: 
Experimental FORC diagram of a Kodak sample. Note the similarity to the FORC 
diagram of the EASG shown in the left diagram. (From I193I V 



6. Equilibrium Properties of Classical XY and Heisenberg 
Spin Glasses 

6.1. Continuous Spin Models in Three Dimensions 

Compared to Ising spin glass models, less effort has been devoted to the 
models with continuous spins. One of the reasons might be that the min- 
imal model of spin glasses, the Ising spin glass model, already shows a 
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highly non-trivial behavior and poses serious difficulties to reach a final 
conclusion. Nonetheless, the continuous spin models are worth studying for 
several reasons. An obvious reason is that, in many real spin glass materials 
such as AuFe and CuMn, the magnetic anisotropy is much smaller than the 
(isotropic) exchange interaction. For these materials, therefore, the Heisen- 
berg spin glass models would be more realistic than the Ising spin glass 
models. 

However, the continuous spin models in three dimensions had long been 
thought of as models without a spin glass transition, which is another rea- 
son why the continuous spin models have not been studied so extensively 
until recently. The scenario that seemed to explain this apparent disagree- 
ment between the models and the materials was that the anisotropy is a 
relevant perturbation and, no matter how small it may be, it changes the 
system to an Ising-like spin glass, giving rise to a ffiiite temperature phase 
transition. This is quite analogous to what happens when a small amount of 
Ising anisotropy was introduced to the isotropic Heiscnberg model (with no 
disorder) in two dimensions. If this scenario is valid, the critical behavior of 
any spin glass material should be the same and fall into the Ising-spin-glass 
universality class. However, there are some qualitative differences between 
experimental results for materials with a small anisotropy and numerical 
results for the Ising spin glasses. In particular in a series of recent torque 
experiments |197ll98j of Heisenberg spin glass samples with varying degree 
of anisotropics indications for a freezing transition even within an external 
field were found (implying the applicability of a mean field or RSB sce- 
nario rather than a droplet scenario to these Heisenberg spin glass systems) . 
More surprisingly, the transition in spin glass materials with lower anistropy 
(more Heisenberg-like) appeared to be more robust against the application 
of an external field than those with higher anisotropy (more Ising like), 
which is in contrast with the traditional view that higher anisotropy implies 
a stronger tendency to order. These observations call for further theoretical 
investigations, and here we report on recent theoretical efforts to scruti- 
nize the existence of a transition in pure Heisenberg spin glasses without 
anisotropy. Below we present two sets of evidences that seem to contradict 
each other. (The conffict has not been solved yet.) While the first set of 
evidences supports a novel phase transition via the chirality scenario (see 
below) , the other suggests an apparently more familiar spin glass transition. 

The first one was proposed by Kawamura and coworkers. (For a review, 
see (199| .') It was suggested that a phase transition may exist in three di- 
mensional Heisenberg model even if there is no anisotropy, and belong to a 
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new universality class different from that of the Ising spin glass model. In a 
series of numerical works (e.g., |71|). the origin of the disagreement between 
the experimental and numerical results was scrutinized by reconsidering the 
role of the magnetic anisotropy in the mechanism of the spin glass transi- 
tion. The resulting hypothesis was that there is a phase transition already 
in the isotropic system but the spins do not show long-range correlations 
even below this critical point. It is therefore difficult, if not impossible, 
to detect an anomalous behavior in spin-spin correlations which had been 
focused on in most preceding numerical simulations. Only the anisotropy 
gives rise to the long range correlations in spins, not the transition itself, by 
coupling them to other degrees of freedom that have been already ordered. 
It was proposed that these degrees of freedom relevant for the transition 
are the chiralities defined as 

Xi = Si ■ {Si+s X Si+s') 

where 5 and 5' are two distinct unit lattice-vectors. 

The dynamical behavior of the three-dimensional Heisenberg model with 
and without a uni-axial spin anisotropy was studica^^ and a clear aging 
phenomena, for both the isotropic and the anisotropic cases, was found in 
the autocorrelation function of the chirality at sufficiently low temperatures. 
For the isotropic Heisenberg spin glass model in three dimensions with the 
Gaussian bond distribution, the critical temperature and the exponent were 
estimated as 

Teh - 0.157(1), /3eh = 1.1(1). 

Such a clear aging phenomena, however, was observed in the bare spin 
degrees of freedom only when the anisotropy exists. It was therefore 
suggestecp2l that this phase transition is not accompanied by a freezing 
of spins in the isotropic case. As can be seen in Fig ll2l the estimate of the 
Binder parameter for the bare spins decreases as the system size increases 
at any temperature and does not show any crossing. In contrast, the Binder 
parameter of the chirality crosses. 

An evidence of a finite temperature phase transition was also found for 
the plane- rotator model (i.e., the XY model with two-component spins) 
in three dimensions. In j76j an equilibrium Monte Carlo simulation was 
performed and a crossing in the Binder parameter defined in terms of the 
chirality was observed while, again, no such crossing was observed for the 
Binder parameter defined with the bare spins. For the 3d plane-rotator 
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Fig. 12. The binder parameter for the ± J Heisenberg spin glass model in three dimen- 
sions. The left panel is for the spin overlap distribution and the right panel is for the 
chirality overlap distribution. (From [721 .') 

model with ± J bond-distribution their estimates are 

Teh = 0.39(3), 1^=1.2(2), 77eh = 0.15(20), Zch = 7.4(10). 

It was further argued that, if a smaU but finite anisotropy exists, the 
spins and the chirahties are coupled, and therefore the spin sector, being 
dragged by anomalies in the chirality sector, also shows anomalies at the 
transition point. It follows that, as a function of the anisotropy, the tran- 
sition temperature continuously approaches its isotropic limit. Moreover, 
the singular part of dragged quantities, i.e., those which diverge when the 
system is anisotropic and do not diverge when it is isotropic, are propor- 
tional to the quantity that drags them, namely, the chirality. The spin glass 
susceptibility, for instance, is given by 

xsG(T,i^)~i?^(r~re)-Tso"' 

The exponent 730 characterizing the spin glass susceptibility is the same as 
the one characterizing the chiral spin glass susceptibility. Their estimates 
of critical indices were compared with those of real materials, and many 
experimental estimates are closer to those of chirality transition than the 
Ising spin glass transition (See, Table 

The second set of evidences contradicts to these findings. In a few recent 
reports on the XY and the Heisenberg spin glasses in three dimensions, it 
was suggested that spins are ordered at the same finite temperature as the 
chiral degrees of freedom even in the isotropic case. The first such evidence 
was obtainecp22l through the computation of the domain- wall energy of the 
XY model in two and three-dimensions. The authors of j2()()| computed the 
spin and the chiral domain-wall energies and estimated the corresponding 
stiffness exponents. In the three-dimensional case, they obtained a positive 
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value for the spin stiffness exponent, 

e^P'" = 0.056(11), 

while they obtained a much larger stiffness exponent for the chiral domain- 
wall energy. The positivity of the stiffness exponent was reconfirmed by an- 
other computation based on the Coulomb gas representationPE] However, 
since the estimated value is close to zero and the actual increase observed 
in the domain-wall energy was only about 10 percent for the whole range 
of system size explored, these results on the stiffness exponent alone could 
not be taken as a conclusive evidence. 

Additional evidences came from Monte Carlo simulations at finite 
temperature p22EnnE2l in which the Binder parameters for the magnetiza- 
tion and the chirality of the three-dimensional ± J Heisenberg model were 
computed. A common crossing of the curves for different system sizes was 
found in both Binder parameters. The two crossing temperatures, one for 
the spin and the other for the chirality, were close to each other. Based on 
this result, the presence of a finite temperature phase transition at which 
both the spin and the chiral degrees of freedom are ordered was suggested. 

This result was confirmed in |74j where again the ± J Heisenberg spin 
glass in three dimensions was studied. There the relaxation of several quan- 
tities at temperatures close to the critical one were measured, starting from 
a completely random initial configuration. In the time regime where the 
time-dependent length scale discussed in subsection 15. II is smaller than the 
system size, the size dependence was absent in the auto-correlation func- 
tion and the observed time-dependence could be regarded as the one in the 
thermodynamic limit. Right at the critical temperature, this time depen- 
dence was well described by an algebraic function for diverging quantities. 
For example, for the spin glass susceptibility, we have 

XSG fx (37) 

with A = 7sG/^i^!^^ Using this asymptotic form, the critical temperature 
as well as the critical exponents (divided by zu) can be determined through 
the observation of the time dependence of some quantity. The exponent ziy 
can be obtained from the temperature dependence of the characteristic time 

t{AT) oc {ATy", 

where AT = T — Tgg. In addition, z was estimated from the results of 
relatively small system sizes for which the correlation length was larger than 
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the system size. In such a setting, the size dependence of the correlation 
time near the critical point should be described by 

T(Ar = 0) cx 

The estimated critical temperature and exponents for the Heisenberg spin 
glass model with ± J bond-distribution in three dimensions were 

Tsg/J-0.21+U-[JJ, 7,g = 1.9(4), 1.1(2), /3 = 0.72(6). 

It should be noted that this estimate of v agrees with some of the experi- 
ments on spin glass materials that appear not to be in the Ising universality 
class. 

A clearer demonstration of the existence of the spin ordering was pre- 
sented in |75j . in which an equilibrium Monte Carlo simulation was per- 
formed covering not only the critical point but also a fairly large region in 
the low temperature phase. Instead of the Binder parameter, which is an 
indirect probe of the presence of the long range order, they used the corre- 
lation length itself as in the estimation of the transition temperature of the 
Ising spin glass model^ (See FigOJ. For both the XY and the Heisenberg 
spin glass with the Gaussian bond distribution in three dimensions, a clear 
crossing in the correlation length divided by the system size as a function of 
the temperature was observed fFig ll3() . The critical indices were estimated 
as 

Tsg = 0.34(2), V = 1.2(2) (XY) 

and 

Tsg = 0.16(2), V = 1.1(2) (Heisenberg). 

It is rather puzzling to have two conflicting sets of numerical evidences, 
represented, e.g., by Figll2l and Fig^J each would pass as a convincing 
evidence if the other did not exist. This puzzling situation has not been 
resolved yet. It is at least clear that there is an unexpectedly large correction 
to scaling and that we need to answer why one of the seemingly legitimate 
methods for detecting the long-range order failed in the present case. 

6.2. Continuous Spin Models in Higher Dimensions 

The Heisenberg spin glass model with Gaussian couplings was studied in 
four and in five dimensions via equilibrium Monte Carlo simulations in 
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Fig. 13. Correlation lengths of bare spins of the XY model (left) and of the Heisenberg 
model (right) divided by the system size. (From 1751 . ) 



As in the study of the three-dimensional cas^^the Binder parameter 
for spins and chirahties was calculated. In the four dimensional case, it 
was found that the spin Binder parameters for various system sizes do not 
intersect while the chirality Binder parameters do. In contrast, in the five 
dimensional case, both kinds of the Binder parameters intersect. Based on 
these results it was suggested that the spin chirality separation occurs only 
in four or lower dimensions while in five or higher dimensions an ordinary 
phase transition, at which both degrees of freedom freeze, takes place. In 
five dimensions, the transition temperature was estimated as = 0.60(2) 
and the critical indices as 

a = -1.0(3), PsG = 0.7(3), 7SG = 1.7(2), vsg = 0.6(2). 

Indications were found that the chiral order phase is characterized by one- 
step replica symmetry breaking unlike the mean-field model. The findings 
for the four dimensional model were not as clear as those for the five di- 
mensional model. In four dimensions, the chiral transition temperature was 
estimated to be Tc = 0.38(2) but no reliable estimates for the critical indices 
for the four dimensional model could be obtained. This may be because four 
dimensions is marginal or nearly marginal and the critical region is rather 
narrow. 

For these numerical evidences based on the Binder parameters, however, 
one has to be aware that the same ambiguities might be present in four and 
five dimensions as in the three-dimensional case discussed in the previous 
subsection. Therefore, while the existence of a phase transition has been 



February 2, 2008 5:33 



WSPC/Trim Size: 9in x 6in for Review Volume 



main 



60 N. Kawashima and H. Rieger 

established in four and five dimensions as well as in three, we still cannot 
make a conclusive statement concerning the nature of the phase transition 
in continuous spin models for the same reason as in the three-dimensional 
case. 

6.3. Potts Spin Glasses 

Potts spin variables are discrete variables that can take on q different states 
Si G {1, 2, .... (7} and that interact with one another in such a way that 
only two cases are discriminated: either both interacting spins are in the 
same state or they are in different states. The Hamiltonian of the g-state 
Potts glass is given by 

^ = - E "^^J- - 1) , (38) 

where the coupling constants Jy are quenched random variables and the 
sum runs over all interacting spin pairs (i.e., all possible spin pairs in a 
mean- field model). 

The g-state infinite range Potts glass with q > displays a dynam- 

ical phase transition at a temperature T^, where the dynamics freezes and, 
for instance, the spin auto-correlations do not decay any more (i.e. they 
reach asymptotically a non- vanishing plateau value). The second, static, 
transition takes place at a lower temperature Tc, below which the EA spin 
glass order parameter has a non- vanishing value and the replica symmetry 
is broken. On the mean-field level there is a close connection between Potts 
glasses for q > 4, and Ising spin glass models with p-spin interactions for 
p > 3, which again display similar self-consistency equations as the mode- 
coupling equations describing structural glasses above the glass transition 
temperature.'^SHEnil Therefore Potts spin glasses can also be regarded as 
simple prototypical models for structural glasses. 

The question unanswered to date is, whether this analytically well- 
established mean-field scenario with the two separate transitions is also 
valid in finite-dimensional models with short range interactions. Recently, 
the g-state Potts glass was thoroughly investigated via Monte Carlo simu- 
lations in three dimensions pSEni on the simple cubic lattice, with q ^ 10 
for a discrete, bimodal bond-distribution as well as a Gaussian distribu- 
tion. In both cases, the first two moments of the distribution were chosen 
such that no ferromagnetic ordering of the Potts spins could occur. It was 
found that for all temperatures investigated the spin glass susceptibility re- 
mained finite, the spin-glass order-parameter remained zero, and that the 
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specific lieat had only a smooth Schottky-like peak. These results could be 
understood quantitatively by considering small but independent clusters of 
spins. These observations imply that there is no static phase transition at 
any nonzero temperature. Consistent with these findings, only very minor 
size effects were observed, which implied that all correlation lengths of the 
models remained very short. 

Moreover, the auto-correlation function C{t) of the Potts spins was com- 
puted. While in the Gaussian model C{t) shows a smooth uniform decay, 
the correlation for the ±J model has several distinct steps. These steps 
correspond to the breaking of bonds in small clusters of ferromagnetically 
coupled spins (dimers, trimers, etc.). The relaxation times follow simple 
Arrhenius laws, with activation energies that are readily interpreted within 
the cluster picture, giving an evidence that the system does not have a dy- 
namic transition at a finite temperature. Hence one can conclude that all 
the transitions known for the mean- field version of the model are completely 
wiped out in three dimensions. 

It should also be mentioned that the mean-field model of the Potts 
spin glass was studied numerically via Monte-Carlo simulations-I^EnnEni 
It turned out the it is extremely difficult to control the finite size effects in 
the mean-field model and that the fourth order cumulant g4{N,T) and the 
Guerra parameter G{N, T) could not be used to locate the static transition 
temperature for the system sizes investigated. Also the spin-autocorrelation 
function C (t) showed strong finite size effects and did not display a plateau 
even for temperatures around the dynamical critical temperature T^. 



7. Weak Disorder 

While there are only very few exactly known facts about spin glass models, 
they serve as helpful guides to numerical studies. This is particularly the 
case in the study of random spin system with weak disorder. An example 
is the findings based on the gauge invariance (a brief review can be found 
in j215| ). A line in the phase space, commonly referred to as the Nishimori 
line (simply N-line, hereafter), was found on which the free energy can 
be computed exactly. The line is associated with the crossover from the 
high-temperature ferromagnetic region to the low-temperature disorder- 
dominated region. The results derived from the gauge invariance seem to 
have a close connection to a number of renormalization group studies. In 
the following, we present a gross survey of these findings and related works. 
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7.1. Phase Diagram of the Discrete Spin Models 

It was pointed oul^^^that the gauge invariance of the partition function un- 
der the transformation Si — > UiSi (with an arbitrary choice ofcr^ G {—1,-1-1} 
for every site) can be used to derive a number of exact relations among var- 
ious quantities in spin glass models. This gauge transformation relates ther- 
mal fluctuations to the geometrical properties of the bond configurations. 
Consequently, a geometric temperature, Tp, can be defined that depends 
only on the bond distribution. This temperature appears in various exact 
relationships derived from this gauge invariance. An example is the identity 
for correlation functions in ± J models: 

[{S,S,)t] - [{S^S,)T{S^S,)TJ (39) 

where {■ ■ ■)t is the thermal average at temperature T. In particular, (• • • 
means the thermal average with the real temperature and the geometrical 
temperature being equal to each other. The symbol [• • • ] denotes, as usual, 
the bond-configuration average. The geometric temperature Tp is defined 
as 

2J 

P={1+ tanh( J/Tp)) /2, or Tp = 

iOgj_p 

where p is the concentration of the ferromagnetic bonds. The higher the 
randomness {p 1/2), the higher the geometric temperature {Tp — > oo). 
The line in the p — T phase diagram defined by Tp = T is the N-line. Along 
this line, (|39|l infers m = q, i.e., the spin-glass order is always accompanied 
by the ferromagnetic order. Therefore, no part of the line lies inside the pure 
spin-glass phase. (Here, we define the spin-glass phase by the conditions 
m = and q > 0.) Another consequence of is that the ferromagnetic 
long-range order is absent at any temperature if the concentration p is such 
that the point {p, Tp) falls in the paramagnetic region in the p—T plane. 
This can be understood by an inequality directly derived from (I39II : 

mS.M < [\{S.S,)t,\] . (40) 

If the point (p, Tp) belongs to the paramagnetic region, the right hand side 
of the inequality is zero in the limit of Rij —f oo, which makes the left 
hand side also zero, meaning the absence of ferromagnetic order at (p, T). 
Assuming that the topology of the T — p phase diagram is as depicted in 
Fig lMf a). this fact implies that any phase in the region p < pn (or Tp > T/v) 
cannot have a ferromagnetic long-range order where pjv is the ferromagnetic 
bond concentration at the N-point. Here the N-point {pn, Tjv) is defined as 
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the point at which the paramagnetic phase boundary intersect with the N- 
hne. It follows that the phase boundary separating the ferromagnetic phase 
from the other low-temperature phase (paramagnetic or spin-glass) must 
be strictly vertical or bend toward the ferromagnetic side. (See Figll4l) 

The inequality H4()|l also means that the ferromagnetic order is maximal 
on the N-line when p is fixed. When Tp < Tjv, as we decrease the temper- 
ature departing from the N-linc, the magnetization decreases and the spin 
glass order parameter becomes larger than the square of the magnetization 
squared. This hints a cross-over from the purely ferromagnetic region to 
the randomness-dominating region at the temperature T — Tp. 

It was further arguec^^ that the phase boundary below the N-line is 
strictly vertical. The argument is based on the fact that the free energy 
on the N-line is the same as the geometric entropy, i.e., the entropy of the 
frustration distribution at p: 

Sip)^-J2PiF)\ogPiF) 

F 

where F is a configuration of frustrated plaquettes and P{F) is the prob- 
ability of F being realized. Since the free energy (or its derivative with 
respect to some external field) must have some singularity at the N-point 
when one moves along the N-line, so does the geometric entropy. The latter 
singularity is, however, solely due to the geometric properties of the bond 
configuration, in which the temperature does not play any role. Therefore, 
it is reasonable to suspect that the same singularity with the geometric 
origin may affect the spin system defined on it at p = regardless of the 
temperature. This leads to a temperature-independent transition point, i.e., 
the strictly vertical phase boundary (See Figll4|l. Note that the argument 
applies also to models with continuous degrees of freedom, such as the XY 
model and the Heisenberg model, if the ferromagnetic phase exists in the 
pure model. It was claimed^"^^ that the vertical phase boundary between 
the ferromagnetic phase and the other is universal. 

In order to substantiate this argument, the authors of |218| considered 
the random Zq model with gauge symmetry, Ti. = — ^{^i — Sj + Jij) 
with Si, Sj, Jij takes on a value 0, 1, 2, • • • or g — 1, and V{- • • ) is a periodic 
fmiction of the period q. This model includes the ± J Ising spin glass model 
as a special case. The sample-to-sample fluctuations of the energy along 
the N-line was computed via the exact relationship between the sample-to- 
sample fluctuations and the thermal fluctuations: 

[{AEf] EE [E'] - [Ef = NbHV'} ~ {Vf) - [{H') ~ {Hf]. 
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Fig. 14. The schematic p — T phase diagram the ± J Ising model for (a) the generic 
case and (b) the two-dimensional case. The thin dashed lines represent the N-line. 

Here, E = {H), Nq the number of nearest-neighbor pairs, and {• • • } de- 
notes a single-bond average defined by {Q{1)} = X^/'Jo 6xp(— 
/ ^^^Q exp(— /3y(Z)). The sample-to-sample fluctuations are proportional 
to the 'geometric specific heat', i.e., the energy diff'erentiated by Kp = l/Tp. 
The results showed that [(Ai?)^] takes its maximum at the N-point. The 
authors of |218j argued that this indicates the existence of a singularity in 
the geometric nature of the system at the N-point, though the nature of 
the singularity is not known possiblly because of the size limitation in their 
calculation. 

It is interesting to note that the correlation length exponent i' at the 
zero-tcmpcrature phase transition from the paramagnetic phase to the fer- 
romagnetic phase in two dimension was founcpi^ED to be very close to (or 
perhaps exactly equal to) the one for two-dimensional percolation. This 
fact, together with the geometric mechanism of the transition, makes it 
very attractive to speculate that the transition is due to a percolation of 
something that is defined geometrically. However, this "something" has not 
yet been identified. 

7.2. Dynamical Properties 

Based on the gauge invariance of the model, several interesting exact rela- 
tions were derivecPsHl fQj- the dynamics. Particularly important is the one 
that relates the auto-correlation function with the totally aligned initial 
state and that with the initialization at T — Tp. Specifically, 



(41) 
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The superscripts specify the initial spin configurations: "(F)" indicates that 
the initial state is the totally aligned ferromagnetic state, and "Tp" an equi- 
librium spin configuration at Tp. The left hand side is the autocorrelation 
function of a spin Si at the t^-th Monte Carlo step (MCS) after the ini- 
tialization (at t = 0) and itself at the {t + tu,)-th step. At i > the system 
evolves with the dynamics of a finite temperature T. The right hand side is 
the autocorrelation function after a quenching from the temperature Tp to 
T (or a sudden heating up, depending upon whether T > Tp or T < Tp). 
Of particular interest is the relation obtained by setting t^, = and T = Tp 
in the above 

This indicates the equivalence of the equilibrium auto-correlation function 
on the N-line to the non-equilibrium relaxation of the single-spin expecta- 
tion value starting from the all-aligned condition. By using this relation, we 
can simply measure the value of a spin without equilibrating at all to obtain 
the equilibrium autocorrelation function. This is a considerable advantage 
from the computational point of view. 

7.3. The Renormalization Group Approach for the Discrete 
Models 

The Harris criteriorP^ is well-known as the criterion by which one can 
decide whether introduction of a weak disorder to a pure system is relevant 
or irrelevant: If the specific heat diverges at the critical point of the pure 
system, i.e., the specific heat exponent a is positive, the universality class 
of the transition in the disordered system will be different from the one of 
the pure system. In particular, it was argued that if the disorder is relevant 
the specific heat divergence is smeared out by the effect of the quenched 
disorder. 

The argument was elaborated further with the help of the renormaliza- 
tion group theory.lSSl The analysis of the randomly diluted m- vector model 
showed that the quenched disorder is relevant around the "pure" fixed point 
(denoted as "P" hereafter), if the specific exponent a is positive at P. As a 
result, the existence of another fixed point, called the random fixed point 
(denoted as "R"), was suggested, and the renormalization group flow leads 
the system to R when starting in the vicinity of P. However, because of the 
small cross-over exponent, i.e., a in this case, the critical region where the 
true critical behavior of R can be observed may be very narrow in many 
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cases, such as the random Ising model in three dimensions (a ~ 0.11). 
For instance results of Monte Carlo simulations of the site-disordered Ising 
systertP^were interpreted in terms of a non- universal behavior, i.e., contin- 
uously varying critical exponent depending on the strength of the disorder. 
It was argued that this apparent non-universal behavior can be interpreted 
as a cross-over from P to R. 

As for the critical indices at R, the following values for three dimen- 
sional systems were reported^Sii based on the four- loop order field-theoretic 
renormalization group calculation: 

P = 0.348, 7 = 1.321, v = 0.671, 77 = 0.032. (42) 

In 12211 the role of the gauge invariance in the renormalization group 
method was scrutinized. It was found that the defining condition of the 
N-line, Tp = T is invariant under a certain renormalization group transfor- 
mation. More specifically, in an extended phase space an invariant manifold 
exists that intersects with the p — T plane with the intersection being the 
N-line. A fixed point is located on this manifold and can be identified with 
the one corresponding to the multi-critical point that is shared by three 
phases: paramagnetic, ferromagnetic and spin-glass. One of the directions 
of the two scaling axes is also found to be parallel to the temperature axis 
while the other is tangent to the N-line. The fixed point is unstable in both 
directions. 

The authors of |226| studied the site-diluted Ising model with Monte 
Carlo simulation. Their estimates for the critical indices were 

P = 0.3546(18), 7 = 1.342(5), ly = 0.6837(24), 77 = 0.0374(36). 

These are very close to the corresponding estimates of the indices for the 
pure model such asP^ 

P = 0.3250(15), 7 = 1.241(2), ly = 0.6300(15). 

Therefore it is technically rather difhcult to discriminate a critical point 
corresponding to the random fixed point from the pure fixed point. 

The numerical renormalization group approach based on Monte Carlo 
simulatiorPSS suggested the simplest flow diagram consistent with all of 
these predictions. The variance and the mean of the domain wall free en- 
ergy were computed and analyzed with the method of the domain-wall 
renormalization groupP^ For the randomly diluted spin systems in three 
dimensions the flow diagram turned out to be the simplest one consistent 
with preceding theoretical predictions (see Fig [TK)l . The RG flow diagram 
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"randomness 



"randomness" 



Fig. 15. A schematic RG flow diagram for (a) non-frustrated models (e.g., the randomly 
diluted ferromagnet) and (b) frustrated models (e.g., the ±J model) in three dimensions. 



for the bond-diluted model consists of two finite temperature fixed points 
(the pure fixed point (P) and the random fixed point (R)) together with 
two zero-temperature fixed points (the ferromagnetic fixed point (F) and 
the zero-temperature percolation fixed point (Z)). The two eigen values of 
the linearized RG transformation were evaluated at R. The positive one 
turned out to be 



yi = 1.47(4) 

which is consistent with the previous estimates of = l/y ' 
above. The other negative one was 



0.68 mentioned 



y2 



-f.3(4). 



A similar calculation for the site-diluted system yielded the result consistent 
with these. 

When in addition the system is randomly frustrated, the simplest flow 
diagram must contain at least three more fixed points: the multi-critical or 
Nishimori fixed point (N), the spin-glass critical fixed point (C), and the 
spin-glass zero-temperature fixed point (SG). The findings in j228j were 
indeed consistent with this simplest flow diagram (Fig ll5|) . In this case the 
exponent at R was estimated as 

2/1 = 1.52(2), and = -0.42(13). 

While the positive one agrees with the estimate for the non-frustrated sys- 
tem, the other does not agree. The reason for this has not been clarified 
yet. 
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7.4. The Location of the Multi- Critical Point 

For the location of the multi-critical point, or the N-point, of the ±J 
model, there are a number of numerical estimates. In particular, the non- 
equilibrium simulation was employed effectively to deduce the equilibrium 
properties. The basic idea of the method can be traced back to the dynam- 
ical scaling hypothesis.'22£-21Q. Various scaling relations, such as l|S7jl . can be 
derived from it. It is assumed that a generating function, which is an ex- 
tension of the equilibrium free energy, exists and that the time-dependence 
of various quantities can be obtained as its derivatives. The hypothetical 
generating function has the form 

/(AT, h, L, t) = L-'^f{ATLy, HLy^tL-') (43) 

analogous to the ordinary form of the finite size scaling. While this form is 
presumably valid in the near-equilibrium time regime, it was suggestecP^ 
that there exists an initial time regime where the dynamics is governed by 
another critical exponent independent of those which appear in (|43(l . 

Using the scaling relations derived from 1)43(1 . the non-equilibrium re- 
laxation of the system was studied. The authors of j231j measured the 
magnetization m{t) = [{Si{t))^''] on the N-line, which is proved to be the 
same as the equilibrium autocorrelation function q{t) as mentioned above 
(see subsection l7.2ll . Then they used the following scaling form 

m{t) = L-^^yq{{T - TN)Ly,tL~') 

that can be derived from ((43(1 . In particular, the long-time asymptotic form 
at T = Tn 

was used to obtain the effective exponent at a finite time, 

_ dlogmjt) 
^^^> ^ ~ dlogt ' 

which should converge to A = /3?//z as t ^ oo at the critical point. A good 
convergence to a finite value was observed around the N-point for the two 
and three dimensional models, although the convergence in the latter case 
seemed more unstable. They quoted the values 

PN = 0.8872(8), A = 0.021(1) (44) 

for the two dimensional model and 



PN = 0.7673(3), A = 0.090(3) 



(45) 
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for the three dimensional model. 

In two dimensions, the exact location of the N-point was predicted by 
|232j . The authors studied the duality transformation formulated in [2'M\\ . 
applying it to a random model with Zg symmetry. The model is defined in 
terms of variables, ^i, each taking one of q values = 0, 1, 2, • • • , q — 1). 
The partition function is defined as 

where V{1) is a periodic function defined on integers with the period q and 
Jij is a quenched random variable that takes on 0, 1,2, •••,(7 — 1. They 
obtained an equation that yields a (possiblly exact) value of the critical 
concentration pc- The resulting values of pc agreed with numerical estimates 
for various models such as the ±J Ising model and the three state Potts 
gauge glass. In addition, a similar equation derived for the random Ising 
model with Gaussian bond distribution yielded a value very close to the 
numerically obtained critical value of J/A J where J and A J are the mean 
and the standard deviation of the bond distribution, respectively. 

The replica method was employed in the derivation of the key equation; 
the duality relation of the n replica system was considered. As reported in 
|234j ■ a quantity, which we denote as a, defined with the physical parame- 
ters, such as p and T for the Ising model, has a very simple transformation 
rule under the duality transformation: 

a q"'/a, 

Therefore, a = q"/^ defines a manifold that is invariant under the duality 
transformation. The virtue of this equation is that we can explicitly take 
its n ^ limit, resulting in 

X:paog('x:e^(''+')-^«) =ilog9, (46) 
;=o \»?=o / 

where pi is the probability of Jy being /. If a well-defined duality equation 
in the n limit exists and if a self-dual point lies on our real p — T plane, 
the condition (|46|l must be satisfied on it. The problem is that we do not 
know what the n — > self-dual equation is nor whether a self-dual point 
lies on our p — T phase diagram. 

However, with the assumption that a self-dual point coincides with the 
N-point, it becomes possible to obtain a number, which may be an exact 
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value of Specifically, the intersection of the N-line and the line defined 
by H46|l for the ± J Ising model is located at 

Pc = 0.889972 • • • 

which agrees with preceding numerical estimates of the N-point such as 

Pc = 0.8905(5) 

in |235| . This remarkable agreement, together with similar agreements for a 
few other models, seems a little too much to regard as a mere coincidence; 
the exact transition points of these models may have been derived from the 
duality. 

7.5. Phase Diagram of the Random XY Model in Two 
Dimensions 

The effect of weak randomness in the XY model in two dimensions is of 
special relevance due to its relationship to various other models, such as 
the Coulomb gas and Josephson junction arrays. The model discussed most 
frequently is the random phase XY model 

= - ^ Jij cos(6'i ~ 9j ~ (t)ij) 
(y) 

with Jij being a uniform constant J > and 0y a quenched random 
variable. However, other models, such as the one with random Jy and 
uniform 0^ j , most likely have essentially the same property. 

In an early studj'^Ssl the random phase model was studied and it was 
suggested that a sufficiently strong disorder destroys the quasi-long-range 
order of the pure XY model. More strikingly, the amount of the disorder 
sufficient to destroy the quasi-order is vanishing as we approach the zero 
temperature. Consequently, the system undergoes a re-entrant phase tran- 
sition as the temperature is lowered with the magnitude of the disorder 
being fixed at a sufficiently small value. 

The latter prediction concerning the re-entrance phase transition was 
corrected later by a number of groupsP '^^ * ^^^ * ^^'^ * ^^'^ * ^^^ * Here we follow the 
heuristic argument given by the authors of |288| . They argued that there 
is a disorder induced phase transition at zero temperature in the two- 
dimensional XY model. Generalizing the Kosterlitz argument for the KT 
transition, they considered the balance between the energy-cost for creating 
a topological excitation (vortex) and the energy-gain due to the interaction 
of the vertex with the disordered potential. They estimated the probability 
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Pdefoct of creating a vortex on a given site being energetically favored. It 
depends algebraically on the area of the system A with an exponent a: 



where a depends on the magnitude of the disorder and the temperature. 
The ordered phase is stable against the introduction of disorder when the 
total number of sites on which the creation of the vortex is favored is zero 
in the thermodynamic limit. That is, Apdcfoct — 0. Therefore, a = 1 defines 
the phase boundary. This result contradicts preceding studies such as the 
one mentioned ahov^^^^ that suggest that the ground state of the pure 
system is unstable at any finite disorder. 

They argued that the preceding renormalization group theories failed 
to capture the correct physics because they neglected the fluctuation in 
the local energy gain due to disorder potential, which can be very large 
with a small but finite probability. In particular, we cannot neglect such 
a fluctuation when it exceeds the thermal fluctuation. This means that 
the previous argument may fail near and below the N-line because the N- 
line can be regarded as the cross-over line below which the geometrical 
fluctuation dominates (see subsection 17. 

They also discussed a finite temperature phase diagram introducing 
the thermal fluctuations in their argument. This yields another interesting 
feature of their results; the straight phase boundary between the ordered 
phase and the disordered phase in the T — Tp phase diagram, in agreement 
with the Nishimori's claim (see subsection I7.1|l . Generalizing the simple 
argument based on the energy-balance, they computed the exponent a for 
finite temperature: 



where T* and A are model-dependent parameters. The phase boundary is 
again given by a = 1. Therefore, for the region where the disorder fluctu- 
ation dominates, {Tp > XT), the phase boundary is Tp — T*, independent 
of the temperature, whereas it depends on the temperature in the region 
where the thermal fluctuation dominates (Tp < AT). If one identifies, quite 
naturally, the line Tp = XT with the N-line, the result is perfectly in parallel 
with Nishimori's picture based on the gauge invariance. 

In fact, the gauge invariance can be used for deriving exact results for 
continuous spin models, such as the XY models and the Heisenberg models. 
In |215j . the author derived a number of exact equations for the Villain 



Pdofcct 



A-' 




(47) 
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model on the N-line. For example, the exact solution for the energy along 
the N-line was obtained. As we have seen in subsection 17. II for the discrete 
models, one can argue, based on the exact results, that the N-line (T — Tp) 
is the cross-over line separating the purely ferromagnetic region and the 
disorder dominant region. Interestingly, the renormalization group theory 
on the XY model^^^ indicates that there is a freezing transition or a cross- 
over at T = 2Tp with the difference of a factor 2 from the N-line. It is 
plausible that this difference is only due to the approximation involved in 
the renormalization group theory, and they both reflect the same physics. 

8. Quantum Spin Glasses 

A quantum spin glass is a magnetic system that can be described by a 
quantum mechanical Hamiltonian with spin-glass like features (randomness 
and frustration). In such a system, a spin glass phase may exist while at 
the same time quantum fluctuations play an important role, possibly a 
dominant role, in particular, in the absence of thermal fluctuations at zero 
temperature. Such a Hamiltonian is, for instance, the spin-1/2 Heisenberg 
spin glass 

(u) 

where tj^'^'^ are Pauli spin-1/2 operators, random exchange interactions 
(e.g., Gaussian), and the sum runs over all nearest neighbors on some d- 
dimensional lattice. Another example is the Ising spin glass in a transverse 
field 

H = -Y,J,,ata^,+TY,af , (49) 

(y) 

where T denotes the transverse field strength. This Hamiltonian becomes 
diagonal if T is zero, in which case it reduces simply to the classical Ising 
spin glass that we have discussed in the previous sections. Thus the role of 
the parameter F is to tune the strength of quantum fluctuations, they do 
not play a role in the equilibrium statistical physics of a diagonal Hamil- 
tonian. An important experimental realization of this model Hamiltonian 
is the system LiHoa:Yi_2;F4,'221 an insulating magnetic material in which 
the magnetic ions (Ho) are in a doublet state due to crystal field splitting. 
The interactions between Ho ions can be described by an Ising model with 
dipolar couplings. For x = 1 the system is a ferromagnet with a critical 
temperature of Tc = 1.53 K at F — and as x is reduced the critical 
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temperature decreases. For concentrations below 25% Ho and above 10% 
Ho a thermal phase transition to a spin glass phase occurs indicated by a 
diverging nonlinear susceptibility (for instance at a: = 0.167 the spin glass 
transition temperature is Tg = 0.13K at F = 0). If a transverse field is 
applied (F > 0) the spin glass transition temperature decreases monotoni- 
cally to zero (see Fig. 116(1 . This particular point, at zero temperature and 
at a critical field strength is what we denote as a quantum-phase-transition 
point.EH 



1.0 



l_ 0,5 



"i^OO O.05 0,10 0.15 

TCK) 

Fig. 16. Phase diagram of LiHoo.i67Yo.833F4 according to the measurement of the 
nonUnear susceptibiUty. From 12431 . 

Earlier reviews on quantum spin glasses and in particular the Ising spin 
glass in a transverse field can be found in |245l246l247j . Here we try to 
focus on a number of new developments that have been made since then. 

8.1. Random Transverse Ising Models 

The generic phase diagram for the EA Ising spin glass model in a transverse 
field F is shown in Fig. 1171 for two dimensions and for three dimensions. In 
the three-dimensional case, starting from the classical spin glass transition 
temperature Tc for F = the critical temperature decreases monotonically 
with increasing transverse field strength F until it reaches T = 0. One 
expects that the universality class of the transition at any non-vanishing 
temperature is the same as the one of the classical Ising spin glass transition 
at Tc- The zero-temperature quantum phase transition, however, establishes 
a new universality class. This transition exists in any dimension, including 
one and two dimensions. A critical value Fc for the transverse field strength 
separates a disordered or paramagnetic phase for F > Fc from an ordered 
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phase for F < r^. This transition is characterized by a diverging length scale 
^ ~ |r — Fcl"'^ and a vanishing characteristic frequency uo ^ AE ~ 
The latter is the quantum analog of "critical slowing-down" in the critical 
dynamics of classical, thermally driven transitions. The new and most im- 
portant property occurring at zero temperature in the random transverse 
Ising model is the infinite randomness fixed point (IRFP) that governs the 
quantum critical behavior at the critical value Tc of the transverse 
One feature of the IRFP is that the dynamical exponent z is formally in- 
finite, the relation between length and energy scales is not algebraic but 
exponential: AE ^ exp{—A^^). 




Fig. 17. Left: Generic phase diagram for the two-dimensional Ising spin glass in a 
transverse field P. Since no spin glass phase is present in d = 2 for T > 0, only a 
quantum spin glass phase and a quantum phase transition at T = exists. Approaching 
the quantum critical point at Fc by decreasing the temperature T, the correlation length 
diverges like T~^/^, where z is the dynamical critical exponent (if z is formally infinite, 
it increases logarithmically). Right: Generic phase diagram of a three-dimensional Ising 
spin glass in a transverse filed. The classical transition temperature (at F = 0) is Tc and 
the corresponding classical correlation length exponent is f^. 



To describe this scenario we generalize the discussion of the transverse 
Ising spin glass by including also random ferromagnetic interactions Jij > 0, 
because many more analytical and numerical results are available for the 
ferromagnetic rather than the spin glass case and the same main features 
are expected to hold in both cases. 

Let us start with a review of the one-dimensional case, in which the 
sign (if it can be negative) of the nearest neighbor couplings can be gauged 
away so that all interactions are ferromagnetic and the resulting model is 
the random Ising chain in a transverse field or a random transverse-field 
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Ising model (RTIM) in one dimension: 



{J^ > 0) 



(50) 



A uniform transverse field is represented hy hi = T for all sites. Since this 
case and the case of random transverse fields turn out to belong to the 
same universality class, we also consider random transverse field here. The 
couplings Ji and the transverse fields hi are random variables with distribu- 
tions 7r( J) and p{h), respectively. The Hamiltonian in H5U|I is closely related 
to the transfer matrix of a classical two-dimensional layered Ising model, 
which was first introduced and studied by McCoy and Wu!^^ Extensive 
researches on this model were initiated by D. Fishei'^^ with an applica- 
tion of the Ma-Dasgupta-Hu renormalization group scheme,^^^ followed by 



numerical and analytical workl ^'^^ l ^-^^ l ^^'^ l ^-^^ l ^^^ l ^^'^ I ^ ^SHEnSl We briefly sum- 



marize the results. The quantum control-parameter of the model is given 



For S < the system is in the ordered phase with a non-vanishing average 
magnetization, whereas the region ^ > corresponds to the disordered 
phase. There is a phase transition in the system at 6 — with rather special 
properties, which differs in several respects from the usual second-order 
phase transitions of pure systems. One of the most striking phenomena is 
that some physical quantities are not self-averaging, which is due to very 
broad, logarithmic probability distributions. As a consequence the typical 
value (which is the value in an frequent event) and the average value of 
such quantities can be drastically different. Thus the critical behavior of 
the system is primarily determined by rare events that give dominating 
contributions to the averaged values of various observables. 

The average bulk magnetization is characterized by an exponent (3, 
which is /3 = 2 — T where r = (1 + \/5)/2 is the golden- mean. The aver- 
age spin-spin correlation function C(r) = [{cr^af^^)]av involves the average 
correlation length ^, which diverges at the critical point as ^ ~ and 
I'av — 2. On the other hand, the typical correlations have a faster decay, 
since ^typ ^ |(5|"'^'='p with z^typ = 1. 

Close to the critical point the relaxation time tr is related to the corre- 
lation length as ~ where z is the dynamical exponent. The random 
transverse- field Ising spin chain is strongly anisotropic at the critical point. 



by 



5 



[\nh]av - [In J] av 



(51) 



var[ln h] -|- var[ln J] 
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(52) 



which corresponds to z = oo. On the other hand the relaxation time is 
related to the inverse of the energy-level spacing at the bottom of the spec- 
trum tr ~ (Ai?)^^. Then, as a consequence of (|52|l . some quantities (such 
as specific heat, bulk and surface susceptibilities, etc.) have an essential 
singularity at the critical point, and the correlation function of the critical 
energy-density has a stretched exponential decay, in contrast to the usual 
power law behavior. 

Away from the critical point in the disordered phase the rare events with 
strong correlations still play an important role, up to the point , 5 — 5g- 
Above this point, all transverse- fields are bigger than the interactions. In 
the region {) < 5 < Sq^ which is called the Griffiths-McCoy phase, the 
magnetization is a singular function of the uniform longitudinal field Hz as 
"^sing \Hz\^/'' , where the dynamical exponent z varies with 5. At the two 
borders of the GrifRths-McCoy phase it behaves as z « ll25x{l + 0{5)^ 
as (5 \ and z = 1 as 5 5q, respectively. 

All these results could be obtained and understood by the application 
of a Ma-Dasgupta-Hu renormalization group scheme, in which strong 
bonds or fields are successively decimated either by elimination of spins 
(in case of large transverse fields) or formation of strongly coupled clusters 
(in case of large ferromagnetic bonds) . With decreasing energy scale A of 
the bonds and fields to be decimated the typical size L of these strongly 
coupled clusters increases as 



defining an exponent ^ that is 1/2 in the random transverse- field Ising 
chain. Such a cluster typically contains 



spins that essentially behave collectively (for instance in response to the 
application of a longitudinal magnetic field H — and thus generating a 
huge contribution to the spin susceptibility). This defines another expo- 
nent (/), which is (1 -f V5)/2 in the RTIM. Finally there is the correlation 
length exponent v that defines the characteristic length scale of spin-spin 
correlations away from the critical point. 

The RG runs into a fixed point that is fully determined by the geometri- 
cal features of the clusters that are generated asymptotically — very much 



llnAji/''' 



(53) 



^^L*'/'(= |lnA|^) 



(54) 
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in reminiscence of the percolation fixed point in conventional percolation. 
This picture is expected to hold also for higher-dimensional RTIMs, and 
even for the spin glass case. Therefore we summarize its essence here. The 
distribution of the random bonds and fields not yet decimated during the 
RG procedure becomes broader and broader. Hence the name, infinite ran- 
domness fixed point (IRFP). It is characterized by the three exponents V', 
and v and the critical behavior of the physical observables is determined by 
them. For instance the correlation function (at criticality) for two spins at 
site i and j with a distance r from each other is simply given by their prob- 
ability to belong to the same cluster of size r: [Cyjav ^ |ri — rj\^^'^'^^'^'^\ 
Other relations follow straightforwardly from this schemd^^ 

lowest energy scale: — In A ~ L"^ 

magnetic moment: /i ~ (— InA)"^ 

average correlations: [Cij]av jr^ — j^^^''"'^'^^ 

typical correlations: — [InC^jJav ~ i^ijl^i ^ 

finite T-susceptibility: x ~ 1^^)20-^/^- 

finite 7f-magnetization: M - (- lniJ)^"*+'^/''' 

Away from the critical point {S ^ 0) the correlation length is finite and its 
average and typical value scale differently: 



average correlation length: ^av ~ 
typical correlation length: ftyp ^ £,lv^ 
spontaneous magnetization: Mq ~ (_^)'^('^-"#"/') 



In spite of the finiteness of the average correlation length away from 
the critical point still arbitrarily large strongly-coupled clusters exist — 
though with an exponentially small probability — leading to algebraically 
decaying correlations in imaginary time. Phenomenologically, one can see 
that as foUows.l^SEnilLet L be the size of a region of strongly coupled spins. 
In a random system in the paramagnetic phase they occur with an expo- 
nentially small probability P{L) oc exp(— AL''). For instance in the diluted 
ferromagnet strongly coupled regions are connected clusters and their prob- 
ability is p^ , where V is the region's volume and p is the site occupation 
probability (0 < p < 1). Then, A is given by A = | lnp| > 0. The spe- 
cial feature of transverse-field Ising systems is that in first order perturba- 
tion theory the gap of a finite system containing L"^ spins is exponentially 
small: Aq ~ exp(— sL"*). An exponentially small gap means an exponen- 
tially large tunneling time, and combining the two observations on cluster 
probability and relaxation time one obtains an algebraical decay for the 



February 2, 2008 5:33 



WSPC/Trim Size: 9in x 6in for Review Volume 



main 



78 N. Kawashima and H. Rieger 

spin-spin correlation function: C(r) = [(cri(r)cri(0))]av ^ t~^/'^ = t^'^I^'^^^ . 
The parameter z(5) = s/dA is called the dynamical exponent in the Grif- 
fiths phase and it varies continuously with the distance from the crit- 
ical point. The consequences, e.g., for the susceptibility are dramatic: 
X{uj = 0) = /g^^dTC(T) cx T-^+d/^{S) which implies that for z > d the 
susceptibility diverges for T —^ even away from the critical point. Since in 
random transverse- field Ising system z{d) grows without bounds for 6—^0 
(and thus merging with the critical dynamical exponent at 6 = 0, which 
is infinite), there is always a region around the critical point, where the 
susceptibility diverges. 

In general the dynamical exponent z{d) introduced above is expected to 
determine all singularities occurring in the Griffiths-McCoy phase close to 
an IRFF248: 

dynamical exponent: z{6) oc 5^'^'^ 

lowest energy scale: A ^ 

finite ff-magnetization: M ^ fj'^/^i^) 

susceptibihty: x(w = 0) - T-i+d/z{S) 

nonhnear susceptibility: Xmi^^ — 0) ^ rp-i+d/3z(S) 
specific heat: c - T'^/^^) 

The last three tables summarize the scaling predictions at and close 
to a IRFP and in Id they have been confirmed many times, analytically 
and numerically:^^^ ^^■^ ^^'^ ^^'"^ •^■'^ In higher dimensions d> 2 the 

randomly diluted Ising-ferromagnet in a transverse field is a show-case for 
a quantum phase transition governed by an IRFP. The site diluted model 
is defined by the Hamiltonian 

H = - e,£j-af a| - F ^ e^af (55) 

and the bond diluted model by 

H^-jY,e,,ata]-rY,af (56) 

where Si and Sij are random variables that take on the values 1 with prob- 
ability p and with probability I — p. Its phase diagram is depicted in Fig. 

CHI 

Along the vertical line starting from the point (p, F) = (pc, 0) up to the 
multi-critical point the transition from the paramagnetic to the ferromag- 
netic phase is a classical percolation transition.^^'^ Denoting the distance 
from the critical point with 5 — Pc~ p the connectivity correlation length 
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Fig. 18. Phase diagram of the diluted Ising ferromagnet in an transverse field F at zero 
temperature T = 0. 

diverges upon approaching the percolation point 1—pc as ~ | (5 1 The 
number of spins M in the percolating cluster at p = pc scales with the linear 
system-size L as M ^ j^DpcT,:^ where Dp^rc is the fractal dimension of the 
percolating clusters. For small values of the transverse field F one expects 
the percolating cluster to be fully magnetized, which implies that the gap 
scales as A ~ exp(— L^i""''=). This means that ip = Dperc in the IRFP sce- 
nario described above. Moreover, the connectivity correlation function at 
the percolation threshold pc decays as C(r) ~ ^-(d-2+r;porc)^ which means 
that the exponent (f> is given by the relation 2(d— ip(j)) — (d—2 + ?7porc)- To 
summarize the exponents characterizing the IRFP in the randomly diluted 
ferromagnet in a transverse field are related to the classical percolation 
exponents (which are exactly known in dimensions d = 2 and in d > 6) via: 

'^ = '^pcrc, Ip^Dpcrc, = (c? + 2 - ?7perc))/-Dporc • (57) 

For the random bond ferromagnet in a transverse field in dimensions 
d > 2 one has to rely only on numerical calculations: In the 2d ferromag- 
netic case quantum Monte-Carlo simulationP^SEsHl provided an evidence 
for an infinite randomness fixed point with ip « 0.42 and (j) ~ 2.1. Later 
a numerical implementation of the Ma-Dasgupta-Hu RG scheme indeed 
provided another evidence for an infinite randomness fixed point^^ with 
= 0.42 ±0.06, = 2.5 ±0.4, which agrees with the QMC estimate within 
the error bars, and v = 1.07±0.15. For random Ising ferromagnets in a 
transverse field the existence of the IRFP dominating the quantum critical 
behavior thus appears to be confirmed for finite dimensions. Strictly speak- 
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ing detailed numerical studies have only be performed for d = 1 and d — 2 
up to now, but there seems to be no strong argument against the existence 
of the IRFP also in higher, finite, dimensions although one expects the nu- 
merical and experimental visibility of the IRFP to diminish for increasing 
dimension d. In the mean field limit {d —^ oo) the quantum phase transi- 
tion is not described by an IRFP and obeys conventional scaling laws. In 
particular z is finite and Griffiths-McCoy singularities are absent. 

What about the spin glass case? Quantum Monte Carlo simulations on 
the Ising spin glass with a transverse-field have been performed for the 
cases d = d^HHEHl and d = 3,^^1123 they are reviewed in EIHEIZI- The 
main result is that the numerical data appeared to be compatible with a 
finite value for the dynamical exponent in d = 2 and 3 and that the crit- 
ical behavior can be described by conventional scaling laws. However, the 
existence of a Griffiths-McCoy phase away from the critical point has been 
uncovered, with a continuously varying dynamical exponent describing the 
singularities of the susceptibility and non-linear susceptibility. In contrast 
to the quantum Monte-Carlo simulations of the random bond ferromag- 
nets no cluster-algorithm could be used in the quantum spin glass case, 
which restricted the system sizes and in particular the temperatures to 
rather small values (note that anisotropic finite size scaling demands that 
the temperature has to decrease exponentially with the system size at a 
quantum critical point described by an IRFP). In addition a homogeneous 
rather than a random transverse field has been used, which causes strong 
cross-over effects and the true asymptotic scaling behavior might be more 
difficult to extract. Therefore it might very well be that the indications 
found for the absence of a IRFP in the 2d and 3d quantum spin glass are 
still pre-asymptotic and that studies using larger system sizes and more 
sophisticated simulation methods could detect evidence for the IRFP also 
here. 

Finally a word about the consequences of the aforementioned theoret- 
ical developments for the experiments. There it was observed that upon 
approaching the quantum critical point the divergence of the non-linear 
susceptibility was drastically suppressed indicating even the absence of a 
divergence at zero temperature. The numerical results, on the other hand, 
hint at a strong divergence of the non-linear susceptibility at the quan- 
tum critical point — even more than the IRFP scenario. Up to now no 
clear reason for the discrepancy has been pinned down. The possibility of 
a second-order transition turning a first-order one at low temperatures has 
been raised,E^ but this possibility can definitely be ruled out for a sys- 
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tern that can be described by the Hamihonian 15U|I that we discussed here. 
We do not think that dipolar interactions of a niagneticaUy diluted sys- 
tem cause substantial modifications of the picture that emerged for short 
range interactions. At this point one cannot rule out the possibility that 
the transverse field Ising Hamiltonian with quenched disorder is simply not 
a sufficiently detailed description of LiHoo.i67Yo.833f?4. 

8.2. Mean-Field Theory 

As a mean-field model of quantum Ising spin glass, we consider the 
Sherrington-Kirkpatrick model in a transverse field 



The first sum is over all pairs of spins and the couplings Jy are quenched 
random variables that obey the Gaussian distribution with zero mean and 
variance J'^/N, where N is the number of spins. T is the strength of the 
transverse field. Although no exact solution has been found for finite T, the 
phase diagram of this model has been well delineated. At zero transverse 
field the transition is the well-known classical transition of the SK model 
at Tc(r = 0) = J. For sufficiently high temperature and/or sufficiently 
large F, thermal and/or quantum fluctuations destroy the spin glass order, 
yielding a paramagnet .12221 For low T and small F one finds a SG ordered 
phase, apparently with broken replica symmetry. ^^"^ Monte Carlo calcula- 
tion, numerical spin summation^^^ and perturbation expansion^^^ in 1/F 
have determined the phase boundary to some precision. As in the classical 
model, the infinite range interactions apparently wipe out the Griffiths sin- 
gularities discussed in the last subsection. The critical behavior along the 
line rc(F) is expected to be the same as the classical critical behavior, i.e., 
the non- linear susceptibility diverges as Xni ^ (T — Tc{T)) ' with 7 = 1, 
the specific heat exponent is a = 0, etc. 

The zero temperature quantum critical point Fc(T — 0) is in a dif- 
ferent universality class and has been studied in |27()I277I278| . The static 
approximation — the approximation usually applied to small field values 
in which the imaginary time correlation function C(r) = ((Ti(r)fTi(O)) is 
assumed to be time independent — is not valid at T = (large fields) and 
the full functional form of C(t) has to be explored. The dynamical self- 
consistency equations obtained via standard manipulationJ2Z212Z2I -v^as ana- 
lyzed at T = at the quantum critical point in |276l277l278j , and it turned 
out that the quantum critical point is located at Fc ~ 0.7 J. At F > Fc (and 




(58) 
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zero temperature) C(r) decays exponentially with r as t ^ c», indicating 
a gap A in the corresponding spectral density; at F = Fc, C(r) decays as 
l/r^, and in the ordered phase, C(t) qea- The Fourier transform of C(t) 
has the form C((jj) ^ const. — \/uj'^ — ISP- for F > Fc, which is responsible for 
the 1 /r^ behavior at Fc and it turned out that the correlation time diverges 
as ~ 1/A ~ [(F - T^y^ ln(F - Fc)]^/^. Thus we can define an exponent 
zv, anticipating anisotropic scaling in space and time in the short range 
model, which takes the value zv = 1/2 in the infinite-range model. Since 
C(t — > oo) — (jEA is the Edwards- Anderson order parameter, we may also 
define ^ea = (J'c — ^iid it was found that /? = 1. At F = Fc one expects 
C(r) ~ r^P/^'^ ^ which is satisfied with the values obtained. The non- linear 
susceptibility diverges as 1/A, which implies with Xni ~ (F — Fc)"''' that 
7 — 1/2. Studying Gaussian fluctuations around the saddle-point solution 
valid for infinite range one find for the correlation length exponent above 
the upper critical dimension (i.e. d > 8) that v — \/A and therefore z = 2. 
Moreover = in mean field theory. The complete collection of critical 
exponents obtained so far in comparison with the classical model (T > 0, 
where we assume to cross the phase boundary under a non- vanishing angle) 
are as follows: 





P 


7 


V 


z 


quantum [T — Q) 


1/2 


1/2 


1/4 


2 


classical (T = 0) 


1 


1 


1/2 





Note that as a consequence of the absence of Griffiths-singularities in mean- 
field models the dynamical exponent z is finite in contrast to the IRFP 
scenario that is supposedly valid for the finite-dimensional models. In a 
longitudinal field one obtains, in analogy to the classical case, an AT man- 
ifold in the T, F, /i phase diagram below which replica symmetry is broken 
and the system is in the SG phase. 

The dynamics of the model (|58|) in the paramagnetic phase has been 
studied in |28()| . where the dynamical single-site self-consistency equations 
have been iteratively solved using a quantum Monte Carlo scheme devel- 
oped in j281j . They mapped the spin-glass transition line in the F-T plane 
using the stability criterion 1 = </xioc; where xioc — 

/o'^ dT C(r) is the 

local susceptibility. They found a second-order transition line ending at 
a quantum critical point at T = in agreement with the argument pre- 
sented above. Going down in temperature to T '--^ O.OIJ and extrapolating 
the results to T = they determined a precise value for the critical field 
Fc = 0.76 ± 0.01, which lies between previous estimates .I^^I^^ The asymp- 
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totic form of C(t) ^ r^^ found in |276| was also confirmed. A comparison 
of the results for the low-frequency susceptibility with the experimental 
curves obtained for LiHoo.i67Yo.833F4 in |243j yields a good agreement. 

A different class of mean-field spin-glass models has been studied in 
|271| — simplified in so far as spherical spins rather than Ising spins were 
considered and more general in so far as p-spin interactions were consid- 
ered. The quantum fluctuations are introduced via a kinetic energy rather 
than the transverse field. The corresponding quantum spherical p-spin-glass 
Hamiltonian is defined by 



where Si are "soft-spins" fulfilling the spherical constraint J2i=i Si{tY — ^ 
for all times t. Quantum mechanics is introduced into the classical p-spin 
glass via the canonical momenta pi that fulfill the commutation relation 
[pi, Sj\ = —ihdij. The multi-spin coupling constants are taken from a Gaus- 
sian distribution with zero mean and variance Jpl/ {2NP~^) with J being a 
constant of 0(1). 

Before we discuss this model we want to clarify the connection to the 
SK model in a transverse field discussed above. The replacement of Ising 
spins Si — ±1 by continuous spins € [— oo, +oo] is often performed 
in the theory of critical phenomena — the discrete symmetry is then en- 
forced by a quartic term sf in the Hamiltonian (this is just the usual 

Ginzburg-Landau theory for critical phenomena with a discrete sym- 
metry) , which also restricts automatically the spin length. Analytically the 
quartic term causes extra complications in all computations, saddle point 
evaluations, RG calculations, dynamical formalism etc. — for which reason 
one often skips it and replaces it by a spherical constraint (either strictly or 
via a Lagrangian parameter having the same effect as a chemical potential). 
Unfortunately the classical spherical mean-field spin-glass model with the 
usual 2-spin interactions does not have a non-trivial spin glass phase. There- 
fore, generalizations to p-spin interactions are sometimes considered.^^ At 
this point a clear connection to the original magnetic system of interest is 
already lost. Nevertheless, one might expect that one can learn something 
about possible scenarios. 

Finally spherical spins cannot be quantized in the same way as Ising 
spins via the application of a transverse field. Therefore they are usually 
quantized via the introduction of a kinetic energy term as in ()6()|l . In addi- 




(60) 
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tion, various analytical techniques available for interacting soft spins with 
kinetic energy, such as the Schwinger-Keldysh formalism, ■^^^ are not avail- 
able for spin operators. The microscopic details of the quantum dynamics 
described by either a transverse field or a kinetic energy term might be very 
different, on large timescales, however, one expects a similar behavior for 
the following reason. To see this, let us consider a model that consists of 
two terms; an arbitrary classical Hamiltonian, Hc\, that is diagonal in the 
z-representation of the spins, and the transverse-field term. Performing a 
Trotter decomposition of the partition function of this model, one obtains 



where L,- is the number of Trotter slices in the imaginary time direction, 
Ar = 13/Lr and K given by e"^^ = tanh(Ar). For At < 1 it is i<r = 
I ln(Arr)|/2. In the last step we neglected a constant factor cosh(Arr)^^. 
If we choose At as a small time cut-off (representing the typical spin fiip 
time) we can approximate the last Trotter sum as the imaginary time path 
integral 



where M = 2KAt = AT|ln(rAT)|. The first term in the integral of the 
action is identical to what one would obtain for the kinetic energy if one 
writes down the imaginary time path integral for the partition sum of the 
Hamiltonian l|6()|l . In this way, the transverse-field term and the kinetic- 
energy term are related. 

In |271| the equilibrium properties of the model were obtained using 
a replicated imaginary-time path integral formalisnpZ2I and analyzing the 
dynamical self-consistency equations for the spin auto-correlation function 
C(t) arising in the limit N ^ oo from a saddle point integration. The result 
for the phase diagram, EA order-parameter and linear susceptibility in the 
case p = 3 are depicted in Fig. El where the parameter F = ti^/{JM) has 
been used — resembling the transverse field strength (since for F ^ one 
recovers the classical case). Above a temperature T* one has a continuous 
transition at a critical point F = ^c{T) from a paramagnetic phase with 
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Fig. 19. Left: Static (thin lines) and dynamic (thick lines) phase diagrams of the p- 
spin model for p = 3. Solid and dashed lines represent second and first order transitions, 
respectively. Right: Magnetic susceptibility (a) and Edwards-Anderson order parameter 
(b) of the p=3 model. (From HTl] ). 

vanishing EA order parameter to a spin glass phase with qea and one- 
step replica-symmetry-breaking (IRSB). Although the EA order-parameter 
jumps discontinuously the transition is second order: there is no latent heat 
(as in the classical case F = 0) and the susceptibility has only a cusp. This 
is due to the fact that the parameter m characterizing the Parisi order 
parameter function q{x) (which is a step function with a single step at 
X — m) is unity at the transition. However, for temperatures below T* 
this parameter jumps at the transition, too, and the transition becomes 
discontinuous; for T < T* the transition is of the first order with latent 
heat and a discontinuous susceptibility (see Fig.lT^. 

8.3. Mean-Field Theory — Dissipative Effects 

An important question that arises for interacting quantum spins at low 
temperatures are the effects of a dissipative environment P'^'^ l ^^ This is 
usually described in terms of its collective excitations, lattice vibrations, 
spin or charge fluctuations, etc., which may be thought of as an ensemble 
of independent quantum harmonic oscillators. A concrete 
example of a single quantum degree of freedom, a spin-1/2 or a so-called 
two- level-system (TLS), coupled to a bath of bosons is the well-known spin- 



boson-modelPam 



H — Hs + Hb + HsB 



(63) 
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where Hs, Hb and Hsb denote the Hamihonian of the system, the bath 
and their coupHng, respectively. These are given by 

Hb ^ -'Y^{pi/m^ + mnUjlxD ^^^^ 



Hsb = - ^ c„x„CT^ 



where T is the transverse field (or tunneling matrix element in the con- 
text of TLSs), n the index enumerating an infinite number of harmonic 
oscillators with coordinates and momenta, Xn and p„, and mass and fre- 
quency, m„ and (jj„, respectively. The constant c„ is the coupling be- 
tween oscillator n and the spin. The spectral density of the environment, 
/(w) = TT ^^(|c„p/(m„a;„)(5(tj — ujn)), is commonly assumed to take the 
standard form^Si. 

/(w) = 2a?i(tj/wph)""^we-"'/'^= , (65) 

where a is a dimensionless coupling constant, lOc a high frequency cut-off 
(which can be set to Wc = oo if < s < 2), and tjph a phonon frequency 
necessary in the non-ohmic (s ^ 1) case to keeps a dimensionless. 

With standard techniqueJ^^^l^ini one can integrate out the oscillator de- 
grees of freedom to express the partition function of the system solely in 
terms of the spin variables 

Z = Tre'^^ = j Va{T)Texp{~S/h) , (66) 

where J 'Da{T) denotes a path integral over all spin configurations (in time), 
T is the imaginary time ordering operator and the action is 

5 = -/ dTTa^{T)-- dr dr' K{t ^ T')cr' {T)a' {t') . (67) 
Jq 2 Jo Jo 

The kernel K{t) is related to the spectral density I{lu) and is for the ohmic 

case (s = 1) essentially an inverse square K{t — r') oc a/{T — r')^. The 

effect of the dissipative environment is therefore a long range interaction of 

the quantum spin in imaginary time. In analogy to the Ising model with 

inverse square interaction^^SJJ depending on the strength of the coupling 

constant a, the system is ferromagnetically ordered or paramagnetic in the 

imaginary time direction; for large a the spin is frozen and for small a the 

spin will tunnel. 

Indeed, for the ohmic case, at zero temperature, there is a phase transi- 
tion at a = ]^ |28b | 287 | Pqj. q, ^ 2 there is tunneling and two distinct regimes 
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develop. If a < 1/2, the system relaxes with damped coherent oscillations; 
in the intermediate region 1/2 < a < 1 the system relaxes incoherently. 
For a > 1 quantum tunneling is suppressed and (tr^) ^ 0, signaling that 
the system remains localized in the state in which it was prepared. These 
results also hold for sub-Ohmic baths while weakly damped oscillations 
persist for super-Ohmic baths.'^^ At finite temperatures (but low enough 
such that thermal activation can be neglected) there is no localization but 
the probability of finding the system in the state it was prepared decreases 
slowly with time for a > ac- 

These conclusions hold for a single spin interacting with a bath. The 
question then arises as to which are the effects of the interplay between 
the spin-spin interactions and the spin-noise coupling in the physics of the 
interacting system. In |282| the effect of a dissipativc bath on a mean- 
field spin glass model with p-spin interactions has been investigated. They 
studied the dissipative spin-boson system H64|l for N interacting spins H — 
Hs + Hb + HsB, where the bath Hamiltonian is the same, the coupling 
Hamiltonian gets an additional sum over the spin index i and Hs is now 
the p-spin Hamiltonian with transverse field 



N 



i=l 



The second term, namely, the multi-spin interaction term is the same as 
the one in (|60(l . For the reason explained in the last section it is analytically 
easier to study spherical spins instead of quantum spin-1/2 degrees of free- 
dom and the quantization of the spherical spins is done via the introduction 
of a kinetic energy term. The partition function then reads 



Z = J Va{T) exp{~S/h) 



(69) 



with the action 

clip 



S 



dT 



i ii<...<ip 

hp M3 



+ zJ2[^Ur)-l] - dT dT' K{t-t')s,{t)s,{t') ,m 
J JO Jo 



where the first term is the kinetic-energy term already motivated in (|61I62I) 
replacing the transverse- field term, the second is the p-spin interaction term, 
the third a term with the Lagrangian multiplier z enforcing the spherical 
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constraint and the last term is the long range interaction imaginary time 
(|67|l that is generated by the integration over the bath variables. 

Starting from 1)70(1 the saddle point equations for the self-consistent 
single-spin dynamics were derivecP^ and the phase diagram computed. 
Analogous to the non-dissipative case discussed in the previous subsection 
a critical line with a second-order section (close to the classical critical point 
{Td,V = 0)) and a first-order section (close to the quantum critical point 
(T = 0, Frf)) was obtained in the presence of a dissipative environment. The 
second order critical line is determined by the condition m — 1, the first 
order critical line is defined as the locus of the points where a marginally 
stable solution first appears with decreasing F for T fixed. For each F and 
a this defines a dynamic transition temperature Td{T,a). The qualitative 
features of the phase diagram, similar to those found for the isolated system, 
see the discussion in the previous section. Notice that the line Td(r, a) lies 
always above T's(F, a), the static critical line that we shall discuss below. 




Fig. 20. Static (left) and dynamic (right) phase diagrams for the p = 3 spin model 
coupled to an Ohmic bath (s = 1). The couplings to the bath are a = 0, 0.25, and 
0.5 from bottom to top. The solid line and line-points represent second and first order 
transitions, respectively. (From 282 ). 

On the right side of Fig. 201 the dynamic phase diagrams obtained for 
p = 3 and three values of the coupling to an Ohmic bath, a = 0, 0.25, 0.5 
is shown. The full line and the line-points represent second and first order 
transition, respectively. 

The first observation that can be made is that in the limit F ^ the 
transition temperature is independent of the strength of the coupling to the 
bath. This is a consequence of the fact that in the limit F ^ the partition 
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function is essentially determined by the zero- frequency components of the 
pseudo-spin which are decoupled from the bath. This result is however non- 
trivial from a dynamical point of view, since it implies that the dynamic 
transition of a classical system coupled to a colored classical bath is not 
modified by the latter. 

The second observation is that the size of the region in the phase space 
where the system is in the ordered state increases with a. Coupling to the 
dissipative environment thus stabilizes this state. This follows from simple 
physical considerations. The interaction term in the action favors spin-glass 
order. Coupling to the bath favors localization and its effect is to reduce 
the effective tunneling frequency. Therefore, in the presence of the bath, 
the value of the bare tunneling frequency needed to destroy the ordered 
state must increase with a. Even if the localized state and the glassy state 
may seem superficially similar, they are indeed very different. In the former, 
the correlation function C(t + t^^t^) approaches a plateau as a function 
of t and never decays toward zero while in the latter the relaxation first 
approaches a plateau but it eventually leaves it to reach zero for t ^ tw 
The fact that the coupling to the environment favors the ordered state also 
reflects itself in the value taken by the order parameters C(t) and qea- As 
a increases, ^^(t) reaches a higher plateau level at long imaginary times. 

8.4. Mean-Field Theory — Dynamics 

The out-of-equilibrium dynamics in real time of the quantum spherical p- 
spin glass model coupled to a dissipative environment, which was discussed 
in the last subsection, was actually studied earlier ■^^^ than the equilibrium 
properties. The response and correlation function are defined in analogy to 
the classical case; C{t + tu,,tw) = N^'^ J2iis'.i't+tw)si{tw) + Si{tw)si{t + tia)) 
(note that the time evolution is now governed by the quantum dynamics 
and C has to be symmetrized in the operators Si{t + t^,) and Si(i^)) and 

i?(t + t^,t^) = iV-l Ss^it + t^)/Sh,{t^). 

In equilibrium the quantum FDT relates R{t) and C(t): 



Away from the critical line, C and R decay to zero very fast with oscil- 
lations. Approaching the critical line Td{a), the decay slows down and if 
> a plateau develops in C. At the critical line the length of the plateau 
tends to infinity. 




(71) 
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In the glassy phase (below the transition) the system does not reach 
equilibrium. For small time differences the dynamics is stationary and time 
translational invariance as well as the QFDT holds: limt^_,oo ^(t+t^, t^) = 
(7+Coq(i). For large times the dynamics is non-stationary, time translational 
invariance nor the QFDT does not hold, and the correlations decay from 
5 to 0. The decay of C becomes monotonic in the aging regime, which 
implies Caging + ^u)) = c{h{t^u) / h{t + tw)) (see subsection 15 .2(1 . One can 
generalize the QFDT in the same spirit as the classical FDT was generalized 
(see subsection 15.411 : 

R{t + t^,,ty,) = 4^{t) / -^e-"^* tanh(X(t + t^, t^)phuj/2)C{t + t^,,uj) 
h J ^vr 

(72) 

with C{t, oj) = 2Re ds exp[iuj{t~ s)]C{t, s). Again, as in the classical case 
(see subsection l5.4(l . Tcff = T/X{t + t^, t^) acts as an effective temperature 
in the system. For a model with two time-sectors it is proposed 



X (^ tyj , t^j ) 



Xst = 1 if t < T{ty,) 
Xagc(?l,T) if t>T{t^) 



with Xago a non-trivial function of h and T and T(tw) is a certain time- 
scale that separates the stationary and aging time-regimes. When t and 
are widely separated, the integration over to in ifT^ is dominated by w ~ 0. 
Therefore, the factor tanh(Xagc(i + i^,, itu)/3fiw/2) can be substituted by 
Xago/3/iw/2 {even ai T = if Xagc(/i,T) = x{h)T when T - 0). Hence, 

Ra.gc{t + twitw) ^ d{t)Xi:,gcl3dt^Ciigc{t + , t^) (73) 

and one recovers, in the aging regime, the classical modified FDT.1221222I 

The self-consistency equations for C{t + tu;,tu,) and R{t + tj^,tuj) were 
evaluated numerically in |292| . An example of the solution is shown in Fig. 
1211 for p = 3. In all figures the following parameters have been chosen: 
zero temperature T = 0, the width of the coupling distribution J = 1, 
the frequency cut-off for the oscillator bath set to Wc = 5, the mass in the 
kinetic energy term M — 1, and the strength of the quantum fluctuations 
H = ah (where a is the spin-bath coupling strength) is ft = 0.1. 

These plots demonstrate the existence of the stationary and aging 
regimes. For t < T{tu]) (e.g. T(40) ~ 5) time translational invariance and 
fluctuation dissipation theorem are established while beyond T{t^) they 
break down. For ft = 0.1 the plateau in C is at q ~ 0.97. C oscillates 
around q but is monotonous when it goes below it. In the inset the depen- 
dence of qEA on ft for T = is presented. Quantum fluctuations generate 
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Fig. 21. Left:The correlation function C(T + tu,,tm) vs r for the p = 3 quantum spher- 
ical p-spin SG model. The waiting times are, from bottom to top, t^ = 2.5, 5, 10, 20, 40. 
Qea ~ 0.97. In the inset, the same curves for t^ = 40 and, from top to bottom, 
h = 0.1, 0.5, 1, 2. Right: The response function for the same model as in the left part. The 
waiting-times increase from top to bottom. In the inset, check of FDT in the stationary 
regime. The full line is , tuj) for t-l-tu, = 40 fixed and £ [0, 40]. The dots are ob- 
tained from Eq.JT^ with Jfst = 1, using the numerical data for Cstat (t) = C(t-|-tm, tm)— 9 
{qea ~ 0.97, see left part). In both cases the response is plotted against t. (From 12921 '). 



a Qea < 1 such that the larger h the smaller qea- The addition of thermal 
fluctuations has a similar effect, the larger T, the smaller qea- In order to 
check the FDT in the stationary regime, in the inset of the right part of Fig. 
1211 a comparison is shown of R{t + t^, t^) from the numerical algorithm for 
t + = 40 fixed and t^j G [0, 40] (fuU fine) with R{t + t^,t^) from Ea. lfT^ 
with X = 1 using Cstat(i) = C{t + tw, tw) — q, q ^ 0.97 obtained from the 
algorithm (dots). The accord is very good if t < T{t^) ~ 5. Finally, when 
one plots parametrically the integrated response x vs. C one finds that for 
C < q ^ 0.97 the X vs C curve approaches a straight line of finite slope 

1/Teff = Xage/T ~ 0.60. 



8.5. Heisenberg Quantum Spin Glasses 

The spin- 1/2 Heisenberg quantum spin glass is defined by the Hamiltonian 
H48|l where the random exchange interactions Jy can be ferromagnetic and 
anti-ferromagnetic. The system cannot be studied efficiently with quantum 
Monte-Carlo methods, due to the sign problem arising from the frustration. 
Therefore, not much is known about these models in finite dimensions, and 
also the mean field theory becomes tractable only in certain limits and 
approximations. 
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8.5.1. Finite Dimensions 

In '294' and later in |295j small clusters of the two-dimensional Heisenberg 
quantum spin glass were studied using exact diagonalization. The average 
total spin in the ground state turned out to scale as S" oc ViV, where N is the 
number of sites. The spin glass order parameter in the ground state extrap- 
olates to a small but non-vanishing value in the thermodynamic limit and 
the spin stiffness does not scale to zero either in the thermodynamic limit. 
Ma-Dasgupta-Hu renormalization group studieJ^SIIllll were performed for 
randomly frustrated spin- ladderpa and in d = 2 and [P^SI for various lat- 
tices and spin- values. The general idea of this RG procedure was already 
described in subsection 18.11 large energies (in the form of exchange inter- 
actions) are successively decimated, ferromagnetic bonds lead to large spin 
formation and anti-ferromagnetic bonds to a spin reduction or even elimina- 
tion in case of equal effective spins connected by the bond to be decimated. 
In two and three dimensions, it was also observed that the final magnetic 
moment to be eliminated increased with system size as ^/N , which cor- 
responds to the aforementioned observation for the ground state spin in 
small clusters. In addition the Ma-Dasgupta-Hu RG calculations in |2H8| 
showed that the probability distribution of the low energy excitations scales 
as P(A) ^ A" with = for 2d and id and that the dynamical critical 
exponent is z = 2 in d = 2 and z = 3/2 in d = 3. 

8.5.2. Mean-Field Model 

The first analytical treatment of the mean-field model of the Heisenberg 
quantum spin glass was performed in j279| applying the replica theory. Al- 
though the solution was confined to the paramagnetic state, the arguments 
for the existence of a low-temperature spin-glass phase were given and the 
critical temperature was estimated. 

Later a Landau theory for quantum rotors on a regular d-dimensional 
lattice was studied in j278j . which is defined by the Hamiltonian 

iJ=|^L2-^J,,n,n,, (74) 

where n.; are M-component vectors of unit length (n^ = 1) and represent 
the orientation of the rotors on the surface of a sphere in Af-dimcnsional 
rotor space. The operators Li^^ {^i < ji.v = 1, . . . , M) are the M [M — 
l)/2 components of the angular momentum L.^ of the rotor: the first term 
in H is the kinetic energy of the rotor with 1/g the moment of inertia. 
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The different components of constitute a complete set of commuting 
observables and the state of the system can be described by a wave function 
^(rii). The action of Li on is given by the usual differential form of 
the angular momentum L^^j, = —i{ni^d / dii, — ni^d/di^). The difference 
between rotors and Heisenberg-Dirac quantum spins is that the components 
of the latter at the same site do not commute, whereas the components of 
the rii do. 

In |278| a Landau theory for this model is derived and it is shown that for 
a suitable distribution of exchange constants this model displays spin- 
glass and quantum paramagnetic phases and a zero-temperature quantum 
phase transition between them. The mean-field phase diagram near the 
T = critical point is mapped out as a function of T, the strength of the 
quantum coupling g and applied fields. The spin glass phase has replica 
symmetry breaking. Moreover, the consequences of fluctuations in finite 
dimensions are considered and above d = 8 the transition turned out to 
be controlled by a Gaussian fixed point with mean-field exponents. Below 
d = 8 a runaway RG flow to strong coupling was found. 

Recently the mean-field Heisenberg quantum spin glass model was gen- 
eralized from the SU{2) spin algebra to an SU{N) symmetry and solved 
in the limit N oo^^ Certain universal critical properties are shown to 
hold to all orders in 1/N. A spin-glass transition was found for all values 
of the spin S and the phase diagram as a function of the spin S and tem- 
perature T was described. The quantum critical regime associated with the 
quantum transition at spin value S ~ and the various regimes in the 
spin-glass phase at high spin are analyzed. The specific heat is shown to 
vanish linearly with temperature. 

The out-of-equilibrium dynamics of the the same model in the same 
limit N —>■ oo, but coupled to a thermal bath, was studied in 300 . It 
was found that the model displays a dynamical phase transition between a 
paramagnetic and a glassy phase. In the latter, the system remains out-of- 
equilibrium and displays an aging phenomenon, which we characterize using 
both analytical and numerical methods. In the aging regime, the quantum 
fluctuation-dissipation relation is violated and replaced over a very long 
time-range by its classical generalization, as in models involving simple 
spin algebras studied previously. 

In the context of Heisenberg spin glasses also the work on metallic spin 
glasses should be mentioned, which were first considered in [30 1| and later 
more extensively in j302j , j303j and [3041 . The main ingredient of a metallic 
spin glass is an itinerant electron systems with random (and frustrated) 
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exchange interactions between the electron spins. Thus in contrast to the 
spin glass systems discussed so far the spins are not fixed to particular sites 
but can diffuse (quantum mechanically) from site to site. These systems are 
motivated by experiments on heavy-fermion compounds such as Yi-^U^Pd 
3 which appear to show a paramagnetic to spin-glass transition with 
increasing doping, x, in a metallic regime. To be concrete the Hamiltonian 
studied in |3U3| is 

where Cia annihilates an electron on site i with spin a =t, i, and the spin 
operator is given by Si^ = X]q/3 '^ta'^ap'^ip/'^^ with the Pauli spin matri- 
ces. The sites i,j lie on a d-dimensional lattice, the hopping matrix elements 
tij are short-ranged and possibly random, and the J^ - are Gaussian random 
exchange interactions, possibly with spin-anisotropies. The remainder _ffint 
includes other possible short-range interactions between the electrons, and 
the resulting total Hamiltonian H has a metallic ground state. 

Starting from this Hamiltonian, in j3D5] , an effective field theory for 
the vicinity of a zero temperature quantum transition between a metallic 
spin glass ("spin density glass") and a metallic quantum paramagnet was 
introduced. Following a mean-field analysis, a perturbative renormalization- 
group study was performed and it was found that critical properties 
are dominated by static disorder-induced fluctuations, and that dynamic 
quantum- mechanical effects are dangerously irrelevant. A Gaussian fixed 
point was found to be stable for a finite range of couplings for spatial di- 
mensionality d > 8, but disorder effects always lead to runaway flows to 
strong coupling for d < 8. Moreover, scaling hypotheses for a static strong- 
coupling critical field theory were proposed. The non-linear susceptibility 
has an anomalously weak singularity at such a critical point. 

In j;-{()4j the competition between the Kondo effect and RKKY inter- 
actions near the zero-temperature quantum critical point of an Ising-like 
metallic spin-glass was studied. In the 'quantum- critical regime,' non- 
analytic corrections to the Fermi liquid behavior were found for the specific 
heat and uniform static susceptibility, while the resistivity and NMR relax- 
ation rate have a non-Fermi liquid dependence on temperature. 

9. Summary and Remaining Problems 

In this review we have tried to provide an overview on the recent develop- 
ments in spin glasses. We concentrated on the topics to which substantial 
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efforts have been devoted in recent years and for which, in our view, the 
most significant progress has been achieved: 1) The numerical investigation 
of the equilibrium thermodynamics of finite-dimensional spin glass models 
with short range interactions using new and powerful methods, such as com- 
binatorial optimization for ground state calculations (through which excited 
states can also be studied) and extended-ensemble Monte Carlo methods 
(e.g., the replica exchange method) for the location and characterization 
of the phase transition. The most challenging task remaining here is still 
the unification of the two paradigmatic pictures: the droplet picture and 
the mean-field scenario. Although being unsolved for several decades now, 
we sec promising steps toward the unification of these two approaches. 2) 
The experimental and theoretical investigation of non-equilibrium dynamics 
and aging phenomena in spin glasses, especially the study of fluctuation- 
dissipation-theorem violations in the glassy phase and the concept of an 
effective (non-equilibrium) temperature. Still a lot of work has to be done, 
experimentally in particular, in order to put the fascinating theoretical 
ideas on a firm and consistent experimental ground. 3) The theoretical 
exploration of quantum effects in spin glasses, statically and dynamically, 
for mean-field models as well as finite dimensional models. Here the most 
demanding challenge appears to be handling the real time quantum dy- 
namics at low temperatures for realistic (i.e. finite dimensional short-range 
interacting) models. 

Unfortunately, we could not cover in this review many "glassy" topics 
related to and inspired by the spin-glass world. First and most actively pur- 
sued in recent years is the theory of the structural glass transition and the 
glass "phase". A lot of progress has been made to make contact between 
spin glasses and structural glasses, the main difference being that the dis- 
order is self-induced in the latter during the freezing process, whereas spin 
glasses live with frozen, time-independent disorder. Recently, it was dis- 
cussed extensively that this might not be a major obstacle in relating the 
two systems. Prom a theoretical point of view, a more substantial difference 
between them with concrete experimental consequences is that spin glasses 
have a well-defined order parameter. A diverging correlation length and 
a divergent susceptibility defined in terms of this order parameter signify 
the transition. In the theory of structural glasses, apparently one cannot 
depend on these helpful vehicles. 

Other related issues are vortex, gauge and Bragg glasses. All three ap- 
pear in the context of disordered superconductors. The gauge glass model is 
essentially an XY spin-glass model with a random (quenched) vector poten- 



February 2, 2008 5:33 



WSPC/Trim Size: 9in x 6in for Review Volume 



main 



96 N. Kawashima and H. Rieger 

tial, originally devoted for describing amorphous granular superconductors 
but later also taken as a paradigmatic model for amorphous high-Tc super- 
conductors. It can be analyzed in the same way as XY spin glasses and also 
has the same order parameter. The major questions concern the existence of 
a finite-temperature phase transition with or without screening effects. The 
vortex-glass model is a model of interacting, elastic magnetic flux lines in 
a random potential, which freezes at low temperatures into a glassy phase 
(which however escapes a clean theoretical description, since an order pa- 
rameter is hard to define). If the disorder is weak, topological defects can be 
neglected and an elastic description is possible, starting from the Abrikosov 
flux line lattice and taking into account its small clastic deformation via 
thermal fluctuations and disorder. For such weakly disordered clastic sys- 
tems a Bragg glass phase is predicted in which the true long-range order of 
the Abrikosov flux line lattice is transformed into a quasi long-range order. 
This glassy phase manifests itself also via an extremely sluggish dynamics. 
Upon increasing the strength of the disorder or the density of the lines (via 
an increased magnetic fleld), topological defects will proliferate, the quasi 
long-range order of the Bragg glass phase vanishes and the system becomes 
a vortex glass. 

Finally one should mention Bose, Fermi and Coulomb glass models, 
which occur in the context of the low-temperature physics of quantum- 
mechanical, disordered, electronic or bosonic systems. The origin of the 
interesting physics is the competition between the quenched random po- 
tential and the interactions between particles, usually long ranged as in 
the Coulomb case. The major issues concern the phase transition between 
conducting (metallic or superconducting) and insulating phases, in which 
the particles are localized. Again the name "glass" for the low-temperature 
phases is justified by the anomalously slow dynamics present here. 

To conclude this tour through the glass zoo we hope to have demon- 
strated that a review on spin glasses should not only be useful for people 
working in the field of frustrated magnets, but also for those encountering 
strong disorder and strong interactions at low temperatures in other fields 
of condensed matter physics. 
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